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-I  This  report  consists  of  three  principal  parts  (Section  2,  3,  and  4);  The 
purpose  of  the  research  reported  in  Section  2  is  to  resolve  the  apparent 
discrepancies  between  different  Physical  Optics  solutions  for  rou^h  sur¬ 
face  scattering.  It  is  also  shown  in  this  section  that  the  appearance 
of  the  so-called  "edge  effect"  is  due  to  premature  truncation  of  a  closed 
surface  integral.  — — 
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v  In  Section  3,  full  wave  solutions  are  developed  for  the  depolarization 
of  the  scattered  radiation  fields  by  conducting  objects  of  Irregular 
shape  above  rough  land  and  sea.  The  research  reported  In  this  section 
deals  primarily  with  the  deterministic  problem.^, 

cIn  Section  4  the  full  wave  approach,  which  accounts  for  both  Bragg  scat¬ 
tering  and  specular  point  scattering  in  a  self-consistent  manner.  Is 
employed  to  evaluate  the  scattering  cross  sections  for  composite  models 
of  rough  surfaces.  It  is  shown  that  the  full  wave  solutions  for  the 
scattering  cross  sections  can  be  expressed  as  a  weighted  sum  of  two  cross 
sections.  The  like  polarized  cross  sections  which  are  associated  with 
the  large  scale  surface  are  independent  of  frequency  and  polarization. 
The  cross  section  which  is  associated  with  the  small  scale  surface  is 
dependent  on  frequency  and  polarization. 
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1.0  INTRODUCTION 


The  recent  impetus  to  produce  rigorous  solutions  to  more  realistic 
models  of  pertinent  propagation  problems  over  a  very  wide  frequency 
range  has  generated  the  need  to  derive  full  wave  solutions  to  problems 
of  radio  wave  propagation  in  dispersive,  inhomogeneous,  anisotropic  and 
dissipative  media  with  irregular  boundaries.  The  considerable  growth  in 
civil  and  military  interest  in  the  development  of  more  reliable  systems 
for  communication,  detection,  navigation  and  positioning,  the  potential 
for  developing  radio  wave  methods  for  remote  sensing  and  the  need  to 
develop  secure  hardened  communication  systems  have  contributed  much  to 
this  renewed  interest.  These  developments  have  been  paralleled  by 
remarkable  advances  that  have  been  made  in  the  availability  of  high 
powered,  very  low  frequency  electromagnetic  sources  as  well  as  the 
availability  of  transmitters  operating  at  optical  frequencies.  The 
ready  access  to  large,  versatile  digital  computers  has  made  it 
possible  to  employ  the  full  wave  approach  to  obtain  numerical  solutions 
to  a  wide  class  of  Important  problems  which  have  hitherto  been  either 
Ignored  or  over-idealized  in  order  to  reduce  them  to  tractable  problems. 

To  perform  the  full  wave  analyses,  it  is  necessary  to  develop 
generalized  field  transforms  that  provide  the  basis  for  the  complete 
expansions  for  the  electromagnetic  fields  in  Irregular  multilayered 
structures  with  varying  thickness  and  electromagnetic  parameters. 

These  complete  expansions  consist  of  the  vertically  and  horizontally 
polarized  radiation  fields,  lateral  waves  and  guided  surface  waves.  The 
generalized  field  transforms  are  used  to  reduce  Maxwell's  equations, 
in  conjunction  with  the  associated  exact  boundary  conditions  for  the 
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electromagnetic  fields,  into  sets  of  first  order  coupled  differential 
equations  for  the  forward  and  backward  traveling  wave  amplitudes. 

The  full  wave  solution^  that  have  been  derived  for  the  scattered 
radiation  fields  from  rough  surfaces  with  arbitrary  slope  and  electro¬ 
magnetic  parameters,  bridge  the  wide  gap  that  exists  between  the 
perturbational  solutions  for  rough  surfaces  with  small  slopes  and  the 
Physical  Optics  solutions. 

Computer  programs  are  being  developed  to  numerically  evaluate 
the  scattering  cross  sections.  Thus  Che  solutions  obtained  using  the 
full  wave  approach  can  also  be  used  by  engineers  who  are  not  necessarily 
familiar  with  the  analytical  techniques  used  in  the  analysis. 

1.1  Stannary  of  Research 

In  his  comprehensive  review  of  theoretical  treatments  of 
scattering  of  electromagnetic  waves  by  rough  surfaces,  Valenzuela  (1968) 
states  that  scattering  models  have  been  developed  "that  make  it  possible 
to  apply  available  theories  to  surfaces  that  cannot  be  treated  purely 
by  perturbation  and  Physical  Optics..."  However,  Valenzuela  notes  that, 
"these  composite  models  are  mostly  based  on  physical  considerations 
and  are  able  to  explain  features  in  radar  cross-section  data  from  the 
sea  that  no  theory  can." 

For  instance,  these  composite  models  sre  used  to  show  that  back- 
scatter  at  near  normal  incidence  is  primarily  due  to  specular  point 
scattering  while  backscattar  at  near  grazing  angles  is  primarily  due 
to  Bragg  scattering.  Thus,  for  perfectly  conducting  rough  surfaces 
the  backscatter  cross  section  for  near  normal  incidence  is  independent 
of  frequency  and  polarization  while  at  near  grazing  angles  it  is 
dependent  on  both  frequency  and  polarization. 


Brown  (1978,  1980)  recently  eaployed  a  combination  of  Burrowa 
perturbation  theory  (1967)  and  Physical  Optica  theory  (Beckmann  and 
Splzzlchino,  1963;  Beckmann,  1968)  to  obtain  a  solution  for  the  back- 
scatter  cross  sections  for  perfectly  conducting  rough  surfaces  in 
terms  of  a  sum  of  two  cross  sections.  The  first  term  in  his  solution 
is  the  specular  point  backscatter  cross  section  associated  with  the 
large  scale  (filtered)  surface  height  and  the  second  term  is  the  Bragg 
scatter  cross  section  associated  with  the  small  scale  surface  height. 
Thus  in  his  work  it  is  necessary  to  decompose  (l.e.,  spectrally  filter) 
the  rough  surface.  To  this  end.  Brown's  specification  of  the  wavenumber 
k^,  where  spectral  splitting  is  assumed  to  occur,  is  based  entirely 
upon  the  characteristics  of  the  small  scale  structure  (1978).  However, 
in  the  approaches  of  Hagfors  (1966)  and  Tyler  (1976),  the  specification 
of  k,  is  assumed  to  be  based  on  the  characteristics  of  the  large  scale 

C 

surface.  Furthermore  Brown's  results  for  the  total  backscatter  cross 
sections  are  critically  dependent  on  the  specific  value  chosen  for  kj 
(Brown,  1978). 

In  this  work  the  full  wave  approach,  which  accounts  for  both  Bragg 
scattering  and  specular  point  scattering  in  a  self-consistent  manner, 
is  employed  to  evaluate  the  scattering  cross  sections  for  composite 
models  of  rough  surfaces  (Bahar,  1981a, b) .  It  is  shown  that  the  full 
wave  solutions  for  the  scattering  cross  sections  can  be  expressed  as 
*  weighted  sum  of  two  cross  sections.  The  like  polarized  cross  sections 
which  are  associated  with  the  large  scale  surface  are  Independent  of 
frequency  and  polarization.  The  cross  section  which  is  associated 
with  the  small  scale  surface  is  dependent  on  frequency  and  polarization. 
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Furthermore  the  specification  of  kg  la  not  baaed  entirely  on  the 
characteristics  of  the  small  scale  structure.  To  determine  the  versatility 
of  the  full  wave  approach,  it  Is  necessary  to  examine  the  sensitivity 
of  the  total  scattering  cross  section  to  variations  in  the  value  of  kg 
where  spectral  splitting  is  assumed  to  occur. 

The  purpose  of  the  research  reported  in  Section  2  is  to  resolve 
the  apparent  discrepancies  between  different  Physical  Optics  solutions 
for  rough  surface  scattering.  It  is  also  shown  in  this  section  that 
the  appearance  of  the  so-called  "edge  effect"  in  Beckmann's  results 
(Beckmann  and  Spizzichino,  1963),  is  due  to  premature  truncation  of  a 
closed  surface  integral. 

In  Section  3,  full  wave  solutions  are  developed  for  the  depolar¬ 
ization  of  the  scattered  radiation  fields  by  conducting  objects  of 
irregular  shape  above  rough  land  and  sea.  The  research  reported  in 
this  section  deals  primarily  with  the  deterministic  problem.  In 
Section  4  these  solutions  are  applied  to  random  rough  surfaces.  It 
is  shown  that  the  full  wave  solution  bridges  the  wide  gap  between  the 
Physical  Optics  approach  and  the  perturbation  solutions. 

Ue  conclude  this  section  with  a  summary  of  the  principal  elements 
of  the  full  wave  approach.  The  principal  properties  of  the  full  wave 
solution  and  its  relationships  to  earlier  solutions  of  scattering 
problems  are  also  summarized  (Bahar,  1981c).  This  summary  is  also 
presented  schematically  in  Figs.  1.1  and  1.2.  The  reader  of  this 
report  who  is  not  familiar  with  the  full  wavs  approach  will  find  this 
summary  useful  even  though  the  details  of  the  full  wavs  method  have 
been  reported  earlier  (Bahar,  1973a ,b,  1974,  1981). 
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TE  SPECTRAL 


Figure  1.1  Principal  Elements  of  Full  Wave  Approach. 
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Figure  1.2'  Principal  properties  of  Full  Wave  Approach. 
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(a)  The  electromagnetic  fields  are  expressed  In  terms  of  complete 
expansions  of  vertically  and  horizontally  polarized  waves.  These  Include 
the  radiation  fields,  the  lateral  waves  and  the  surface  waves  (Bahar, 
1973a, b;  1974). 

(b)  Exact  boundary  conditions  are  imposed  at  the  Irregular  surface. 

(c)  Using  the  orthogonal  properties  of  the  basis  functions  appearing 
in  the  complete  expansions  of  the  fields  .Maxwell's  equations  are 
integrated  over  the  transverse  plane  (y,z)  (Bahar,  1973a, b;  1974). 

Green's  theorems  are  used  to  avoid  term-by-term  differentiation  of  the 
field  expansions. 

(d)  Maxwell's  equations  for  the  electromagnetic  fields  are  con¬ 
verted  into  coupled  first  order  ordinary  differential  equations  for 
the  forward  and  backward  traveling  wave  amplitudes  which  are  only 
functions  of  the  variable  x  (Bahar,  1973a, b;  1974).  (  In  view  of  the 
integration  in  the  transverse  plane  (y,z)  the  telegraphists'  equations 
are  only  functions  of  x) .  The  coupled  equations  for  the  wave  amplitudes 
are  referred  to  the  generalized  telegraphists'  equations  (Bahar,  1981c). 

(e)  Second  order  iterative  solutions  for  the  radiation  fields  are 
obtained  from  the  telegraphists'  equations  on  neglecting  multiple 
scattering  from  the  rough  surface.* 

(f)  A  variable  coordinate  system  that  conforms  with  the  local 
features  of  the  Irregular  boundary  is  introduced  and  the  resulting 
solutions  for  the  scattered  fields  are  shown  to  be  invariant  to 
coordinate  transformations. 

(g)  The  full  wave  solutions  are  also  shown  to  satisfy  the  reciprocity 
realizability  and  duality  relationships  in  electromagnetic  theory. 

*The8e  second  order  iterative  solutions  account  for  wave  scattering  in 


arbitrary  directions. 


(h)  The  full  wave  approach  not  only  accounts  for  scattering  and 
depolarization  of  the  radiation  fields  but  also  accounts  for  coupling 
between  the  surface  waves  ,  the  lateral  waves  and  the  radiation  fields. 

The  versatility  of  the  full  wave  approach  is  demonstrated  by 
determining  its  relationship  to  earlier  solutions.  Thus  on  using  a 
stationary  phase  approach  to  evaluate  the  integrals  for  the  scattered 
fields,  the  full  wave  approach  is  shown  to  reduce  to  the  geometric 
optics  solutions  (Bahar,  1981c) .  On  the  other  hand,  if  the  vector  ii 
normal  to  the  rough  surface  is  replaced  by  its  value  at  the  specular 
points  ns>  the  full  wave  expressions  for  the  scattered  fields  are  shown 
to  reduce  to  the  Physical  Optics  solutions.  Thus  the  Physical  Optics 
approach  is  valid  only  if  the  contributions  to  the  scattered  fields 
come  primarily  from  specular  points  on  the  rough  surface.  Details 
of  this  aspect  of  the  research  are  given  in  Section  3.6. 

In  a  survey  of  the  technical  literature  one  finds  several  different 
forms  of  Physical  Optics  solutions.  The  discrepancies  between  the 
different  Physical  Optics  solutions  and  the  appearance  of  the  so-called 
"edge  effect"  have  been  shown  to  be  the  result  of  premature  truncation 
of  the  closed  surface  integrals.  Details  of  this  research  are  reported 
in  Section  2. 

If  one  assumes  that  the  scale  and  the  slopes  of  the  rough  surface 
are  small,  it  is  shown  that  the  full  wave  solutions  reduce  to  the 
perturbation  solutions.  Details  of  this  aspect  of  the  research  are 
given  in  Section  3. 

The  Physical  Optics  solutions  for  the  backscattered  fields  become 
singular  for  near  grazing  angles.  Thus  in  this  case,  even  if  the  rough 
surface  satisfies  the  radii  of  curvature  criteria  (associated  with 
the  Kirchhoff  approximations  of  the  surface  fields),  the  Physical  Optics 
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solutions  cannot  be  used.  This  is  because  of  near  grazing  angles, 
the  principal  contributions  to  the  backscattered  fields  do  not  come 
from  specular  points  of  the  rough  surface.  (In  this  case  specular 
points,  if  they  existed,  would  be  on  vertical  portions  of  the  rough 
surface).  It  is  shown  that  the  full  wave  solutions  for  the  back- 
scattered  fields  remain  valid  as  one  approaches  grazing  angles 
(Bahar,  1982). 


1.2  Interim  Technical  Reports 


Preprints  of  the  following  manuscripts  were  submitted  to  the 
Contract  Monitor  and  published  as  a  U.S.  Air  Force  Interim  Technical 
Report,  RADC-TR-81-204 ,  August  1981. 

Full  Vave  Analysis  for  Scattering  Cross  Sections 

Part  1  Random  Rough  Surfaces 
Part  2  Composite  Surfaces 


1 ■ 3  References 

1.  Bahar,  E. ,  "Depolarization  of  Electromagnetic  Waves  Excited  by 
Distribution  of  Electric  and  Magnetic  Sources  in  Inhomogeneous 
Multilayered  Structures  of  Arbitrarily  Varying  Thickness-Generalized 
Field  Transforms,"  Journal  of  Mathematical  Physics,  14(11) ,  pp.  1502- 
1509,  1973. 

2.  Bahar,  E. ,  "Depolarization  of  Electromagnetic  Waves  Excited  by 
Distributions  of  Electric  and  Magnetic  Sources  in  Inhomogeneous 
Multilayered  Structures  of  Arbitrarily  Varying  Thickness-Full  Wave 
Solutions,"  Journal  of  Mathematical  Physics,  14(11) ,  pp.  1510-1515, 
1973. 


10 


3.  Bahar,  E. ,  "Depolarization  In  Monuniform  Multilayer ad  Structures — 
Full  Wave  Solutions,"  Journal  of  Mathematical  Physics,  15(2) , 

pp.  202-208,  1974. 

4.  Bahar,  E. ,  "Scattering  Cross  Section  from  Random  Rough  Surfaces — 

Full  Wave  Analysis,"  Radio  Science,  16(3) ,  pp.  331-341,  1981a. 

5.  Bahar,  E.,  "Scattering  Cross  Sections  for  Composite  Random  Surfaces — 
Full  Wave  Analysis,"  Radio  Science,  16(6) ,  pp.  1327-1335,  1981b. 

6.  Bahar,  E. ,  "Full  Wave  Solutions  for  the  Depolarization  of  the 
Scattered  Radiation  Fields  by  Rough  Surfaces  of  Arbitrary  Slope," 

IEEE  Transactions  on  Antennas  and  Propagation,  AP-29(3) ,  pp.  443- 
454,  1981. 

7.  Bahar,  E. ,  "Scattering  and  Depolarization  by  Rough  Surfaces  Near 
Grazing  Angles — Full  Wave  Solutions,"  IEEE  Transactions  on  Antennas 
and  Propagation,  AP-30(4) ,  pp.  712-719,  1982. 

8.  Beckmann,  P. ,  The  Depolarization  of  Electromagnetic  Waves,  Golem 
Press,  Boulder,  Colorado,  1968. 

9.  Beckmann,  P. ,  and  A.  Spizzichino,  The  Scattering  of  Electromagnetic 
Waves  from  Rough  Surfaces,  MacMillan,  New  York,  1963. 

10.  Brown,  G.  S. ,  "Backscattering  from  a  Gaussian-Distributed  Perfectly 
Conducting  Rough  Surface,"  IEEE  Transactions  on  Antennas  and 
Propagation,  AP-26(3) ,  pp.  472-482,  1978. 

11.  Brown,  G.  S. ,  "Correction  to  Backscattering  from  a  Gaussian- 
Distributed  Perfectly  Conducting  Rough  Surface,"  IEEE  Transactions 
on  Antennas  and  Propagation,  AP-28(6) ,  pp.  943-946,  1980. 

12.  Burrows,  M.  L. ,  "On  the  Composite  Model  for  Rough  Surface  Scattering," 
IEEE  Transactions  on  Antennas  and  Propagation,  AP-21 (2) .  pp.  241-243. 
1967. 


11 

13.  Hagfors,  T. ,  "Relationship  of  Geometric  Optica  and  Autocorrelation 
Approach  to  the  Analysis  of  Lunar  and  Planetary  Radar,"  Journal 

of  Geophysical  Research,  71,  pp.  379-383,  1966. 

14.  Tyler,  G.  L. ,  "Wavelength  Dependence  in  Radio-Wave  Scattering  and 
Specular  Point  Theory,"  Radio  Science,  11(1),  pp.  83-91,  1976. 

15.  Valenzuela,  G.  R. ,  "Scattering  of  Electromagnetic  Waves  from  a  Tilted 
Slightly  Rough  Surface,”  Radio  Science,  3/11),  pp.  1051-1066,  1968. 


12 


2.0  RESOLUTION  OF  THE  DISCREPANCIES  BETWEEN  DIFFERENT  PHYSICAL  OPTICS 

SOLUTIONS  FOR  ROUGH  SURFACE  SCATTERING 

2.1  Background 

The  purpose  of  the  analysis  presented  in  this  section  is  to  resolve 
the  discrepancies  between  the  different  high  frequency!  Physical  Optics 
expressions  for  the  scattered  field  derived  by  several  researchers  in 
the  field.  The  analysis  also  resolves  the  questions  associated  with  the 
"edge  effect"  that  appears  in  some  of  the  earlier  solutions.  It  is 
shown  that  the  Physical  Optics  approach  is  only  applicable  for  specular 
point  scattering  and  that  the  so-called  "edge  effect"  which  appears  in 
earlier  derivations  of  the  Physical  Optics  solution  for  rough  surface 
scattering  is  a  result  of  premature  truncation  of  the  closed  surface 
integral  expression  for  the  scattered  fields.  Therefore,  this  term  must 
be  suppressed  even  when  it  is  not  very  small  compared  to  the  scattered 
field  in  the  off  specular  direction.  Since  the  Kirchhoff  approximations 
for  the  surface  fields  are  used  in  the  Physical  Optics  approach,  it 
cannot  account  for  wave  diffraction  by  edges. 

The  Physical  Optics  solution  derived  here  for  arbitrary  source 
excitation  is  shown  to  satisfy  reciprocity  and  realizability  relation¬ 
ships  in  electromagnetic  theory.  The  Integrand  in  the  integral  expression 
for  the  scattered  field  is  identified  with  the  specific  reflectance  (per 
unit  area)  of  the  rough  surface.  Although  the  scalar  acoustic  problem 
is  considered  here  in  detail,  the  results  are  also  applicable  to 
electromagnetic  scattering. 
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The  high  frequency  approximation  of  the  full  wave  solution  for 
the  scattered  field  is  in  agreement  with  the  Physical  Optics  solution 
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derived  here. 

2.2  Discussion 

Applying  Green's  theorem  to  problems  of  rough  surface  scattering 
and  using  the  Xlrchhoff  approximations  for  the  fields  on  the  surface. 
Physical  Optics  (high  frequency)  solutions  for  the  scattered  fields  have 
been  derived  in  terms  of  an  integral  over  the  rough  surface  (Beckmann  and 
Spizzichino,  1963) .  On  performing  integration  by  parts,  a  boundary  term 
which  does  not  depend  upon  the  shape  of  the  rough  surface  appears  in  the 
solution.  When  the  scattering  surface  is  large  (in  terms  of  wavelengths), 
this  term  which  is  identified  by  Beckmann  as  an  "edge  effect"  is  considered 
to  be  very  small  compared  to  the  remaining  surface  integral  and  therefore 
negligible. 

Several  researchers  in  the  field  have  clearly  demonstrated  however, 
that  in  general  the  so-called  "edge  effect"  is  not  necessarily  very 
small  in  comparison  with  the  remaining  surface  integral.  It  has  been 
suggested  that  the  "edge  effect"  ....  "may  become  important  when 
considering  scattering  in  directions  off  specular"  (Fung,  Moore  and 
Parkins, 1965) .  Wetzel  (1966)  states  that  Beckmann's  argument  that  the 
integrated  terms  (edge  effect)  can  be  discarded  "leads  to  an  incorrect 
result  for  the  flat  surface  limit  and  suggests  a  scheme  to  rectify  the 
discrepancy  by  retaining  the  integrated  terms."  Furthermore,  Valenzuela 
(1968)  finds  "that  Beckmann's  approximation  in  which  he  neglected  the 
'Edge  effect'  leads  to  incorrect  predictions  at  angles  away  from  normal 
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incidence"  end  derives  solutions  "without  neglecting  the  'edge  effect' 
terns."  However,  it  is  shorn  (Section  2.4)  thet  when  the  tern  identified 
with  the  "edge  effect"  is  included,  the  Physicel  Optics  solution  (Beck- 
mann  end  Spizzichino,  1963)  does  not  setisfy  the  reciprocity  relationship 
in  electronagnetic  theory. 

In  addition,  a  peculiar  non-vanishing  term  appears  in  the  expression 
for  the  scattered  field  even  as  the  reflection  coefficient  for  the  rough 
surface  is  set  equal  to  zero.  This  unacceptable  result  becomes  more 
pronounced  in  the  off  specular  direction. 

To  compound  the  problems  one  encounters  with  the  so-called  "edge 
effect,"  Wetzel  (1975)  points  out  that  the  three  formulations  for  the 
scattered  field  derived  by  Eckhart  (1953),  Parkins  (1967),  Boyd-Deaven- 
port  (1973)  "give  three  different  results  for  the  scattered  field."  What 
is  not  so  obvious  (according  to  Wetzel)  is  "why  this  should  be."  Wetzel 
goes  on  to  suggest  that  "These  inconsistent  results,  .  .  .  must  somehow 
be  associated  with  the  use  of  the  Kirchhoff  approximation  in  the  scalar 
Helemholtz  integral."  All  three  formulations  cited  above  are  based  on 
the  Kirchhoff  approximations  and  the  discrepancies  are  obvious  even 
before  the  integration  by  parts  and  the  appearance  of  the  term  associated 
with  the  "edge  effect." 

Several  questions  therefore  need  to  be  answered.  Should  the  so- 
called  "edge  effect"  be  included  in  the  solution  if  it  is  not  very  small 
compared  to  the  scattered  field?  If  the  "edge  effect"  does  "become 
important"  and  is  included  in  the  solution,  how  can  one  account  for  the 
fact  that  the  resulting  Physical  Optics  solution  does  not  satisfy 
reciprocity  and  does  not  vanish  as  the  reflection  coefficient  becomes 
vanishingly  small?  How  can  one  resolve  the  discrepancies  in  the  formu¬ 
lations  by  Eckart  (1953),  Parkins  (1967)  and  Boyd-Deavenport  (1973) 
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even  though  all  apply  Kirchhoff  approximations  and  none  Involve  the  so- 
called  "edge  effect"? 

To  this  end,  in  Section  2.3,  the  general  form  of  the  Physical  Optics 
solution  for  the  scattered  fields  is  presented  in  terms  of  a  closed  sur¬ 
face  integral  and  an  arbitrary  source  distribution  (spherical  waves 
rather  than  plane  wave  excitation)  is  assumed.  In  Section  2.4  the  limit¬ 
ing  forms  of  the  Physical  Optics  solutions  with  and  without  the  "edge 
effect"  are  derived  for  plane  wave  excitations.  Realizability  and  reci¬ 
procity  relationships  in  electromagnetic  theory  are  examined  and  the  limiting 
case  in  which  the  reflection  coefficient  vanishes  is  also  considered. 

In  Section  2.5  it  is  shown  that  the  discrepancies  between  the  different 
formulations  for  the  scattered  field  and  the  problems  that  arise  with  the 
so-called  "edge  effect"  in  the  Physical  Optics  solution  is  due  to  the  pre¬ 
mature  truncation  of  the  closed  surface  integral  encountered  in  the  appli¬ 
cation  of  Green' 8  theorem.  On  applying  the  divergence  theorem  (in  two 
dimensions),  the  term  associated  with  the  "edge  effect"  is  shown  to  vanish 
identically  for  all  scatter  directions.  The  integrand  in  the  final  expres¬ 
sion  derived  here  for  the  Physical  Optics  scattered  field  is  identified 
with  the  specific  reflectance  (per  unit  area)  of  the  rough  surface. 
Furthermore,  the  integrand  is  not  proportional  to  the  gradient  of  the 
tough  surface  and  is  therefore  not  undefined  at  edges.  The  Integration 
can  be  performed  only  over  those  portions  of  the  rough  surface  with  a 
non-vanishing  reflection  coefficient.  Shadowing  effects  are  considered 
and  the  distinction  between  the  scattered  field,  the  reflected  field  and 
the  shadow  forming  field  are  taken  into  account  (Morse  and  Feshbach , 

1953). 


16 


Th«  final  form  for  the  Physical  Optics  scattered  field  satisfies 
the  reciprocity  and  realizability  relationships.  It  vanishes  as  the 
reflection  coefficient  vanishes  at  the  rough  surface  for  all  incident 
and  scatter  angles.  Using  the  Physical  Optics  approach,  the  term  asso¬ 
ciated  with  the  "edge  effect"  must  be  suppressed  even  for  scattering 
far  off  the  specular  direction.  The  Phyaical  Optics  solution  which  employs 
a  Kirchhoff  approximation  for  the  surface  fields  cannot  account  for  wave 
diffraction  by  edges. 

In  Appendix  (2. A),  it  is  shown  that  on  following  the  analytical 
procedures  of  Section  2.5,  the  Eclcart ,  Parkins  and  Boyd-Deavenport  solu¬ 
tions  are  in  agreement  provided  the  integration  is  over  the  closed  sur¬ 
face  Aq.  Furthermore,  the  integrand  in  all  the  above  Physical  Optics 
solutions  (which  employ  the  Kirchhoff  approximations)  are  shown  to  be 
equal  at  the  specular  (stationary  phase)  points. 

While  the  "Helmholtz  Integral  Equation"  (2.10)  (associated  with 
the  scalar  problem)  is  considered  here  ir.  detail,  it  is  interesting  to 
note  that  similar  difficulties  arise  when  the  "Stratton-Chu  Integral 
Equation"  (associated  with  the  corresponding  vector  problem)  is  applied 
to  surfaces  chat  are  not  closed  (Barrick,  1965) . 

2.3  Formulation  of  the  Problem 

Consider  the  solution  to  the  scalar  acoustic  potential  ifi(r)  that 
satisfies  the  inhomogeneous  wave  (Helmholtz)  equation 

(V2+  k2)*(?)  -  -Air  s(r)  (2.1) 

in  which  s(r)  is  proportional  to  the  acoustic  source  density  and  the 
acoustic  wave  number  is 

k  ■  w/c  ,  (2.2) 


in  which  c,  the  acoustic  velocity,  is  related  to  the  medium  density  p 
and  elasticity  X  and  acoustic  Impedance  Z  through  the  equations 
(Brekhovskikh, 1960) 
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c  -  ^  ,  Z  -  pc  -  /pi  (2.3) 

The  acoustic  pressure  P  and  particle  velocity  V  are 

V  -  -grad  \|i  ,  P  «■  -iup  \p  (2.4) 

and  the  boundary  conditions  at  an  Interface  between  two  media  charac¬ 
terized  by  c,p  and  c^,p^  respectively  are 

P(?s)  -  P1(rs)  and  n'V(rg)  -  n*V1(rg)  (2.5) 

in  which  n  is  the  unit  vector  normal  to  the  interface.  Thus  at  r  ■  r 

s 

p^  -  P14*1  (2.6a) 


and 


‘  oi  -  ur  M  3*1 

n-V*  -  n-V^  or  ^ 


(2.6b) 


To  facilitate  the  solution  to  (2.1)  consider  also  the  Green's  func¬ 
tion  G(r|r  )  that  satisfies  the  equation 
o 

(V2  +  k2)G(r|r  )  -  -An6(r-r  ) 
o  o 


(2.7) 


in  which  6(r-r  )  is  the  Dirac  delta  function.  For  an  infinite  medium 
o 


characterized  by  parameters  p,c  (Morse  and  Feshbach,1953) 

1 


G(r |rQ)  -  G(Fo|?)  -  eik|r"rol 


r-r 


(2.8) 


in  which  an  exp(-iwt)  time  dependence  is  assumed.  Applying  Green's 
second  theorem  to  the  volume  Vo  bounded  by  the  closed  surface  Aq 
separating  the  two  media,  p,c  and  ,c^  one  gets 


j  0J/V2G-CV2iJ0dVo 


G-GV  ^]*dAo 


(2.9) 


I 
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Thus  (Mora*  and  Feshbach,1953) , 


*<?>  -  |«(to)G(r|to)<lVo  +  ji  - 

-c(f|F;)7o*(p)>;  a*, 


(2.10) 


In  which  both  tha  source  end  observation  point  are  la  the  volume  end 

■JX  -  -ndA  where  n  la  tha  unit  nornel  to  tha  surface  A  pointing  towards 
o  o  o 

the  madiua  in  which  the  source  and  observation  points  are  located  (see 
Fig.  2.1).  If  mediua  p,c  is  infinite  (p-p^  c-c^)  the  surface  integral 
vanishes  and  4*(r)  reduces  to  the  incident  unperturbed  field  \|>*(r) 


(Morse  and  Feehbach>1953).  Thus 


-  velkl?-?ol 


i  -  f  -  ax*'r 
^(r)  -  «<ro>— r 

*  r-r 


If  the  source  s(r  )  is  in  the  vicinity  of  the  origin 
o 

|r-ro|  “  ri-ro»af  ,  ar  -  r/r  -  n1  ,  |r|  -  r1 


the  far  field  expression  for  tp*  is 

.  ikr*  r 

V( r)  <*  ■*—£-  Js(r  )e 
r  ' 


-ikr  *n*  .  ikr* 

•  °  dv0  •  *<“*>  S^T 

r 


(2.11) 


(2.12) 


(2.13) 


in  which  the  gain  function  g(n  )  depends  upon  the  direction  of  the  unit 
vector  n1(6,^)  pointing  from  the  source  (in  the  vicinity  of  the  origin) 
to  the  point  at  r  (Jordan  and  Balsa in, 1968) .  Since  g(n*)  is  not  a  function 
of  r*,  the  distance  fros  the  source  to  the  field  point,  tha  far  field 
approximation  for  the  incident  velocity  vector  5^(r),  (2.4), is  given  by 

V1  -  -V^Sr)  -  -(Ik  -  -i-)nV(r)  *  -ik»V(r)  (2.14a) 


Thus  for  the  far  field 


V1  -  ik  n1  FViup  -  n^/Z 


(2.14b) 


in  which  Z  is  the  acoustic  impedance  (2.3).  Expressing  the  total  acoustic 
potential  <p  as  the  sum  of  the  incident  (unperturbed)  and  scattered 


~V- v  ■ 


mmrnm 
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potential  ip1  and  ifi8  respectively 


tp  -  ip1  +  ip* 


(2.15) 


equation  (2.10)  reduces  to  (Morse  and  Feshbach  1953) 

*8<?)  "  4?  '  G(?l?o)V*<?o)],“dAo 


(2.16a) 


and 


(2.16b) 


The  Kirchhoff  (Physical  Optics)  approximations  for  (|/(r8)  and  VQ^(r®)  are 


*(r8)  -  (1  +  R)4>*(r8)  ,  VQ^(r^)  -  ito^l  -  R)^8) 


(2.17) 


-i  -» 


in  which  R(-n  *n)  is  the  acoustic  reflection  coefficient  at  a  flat  inter¬ 
face  between  medium  p,c  (where  the  incident  wave  originates)  and  medium 
(Brekhovskikh,1960;  Ishimaru ,1978) .  Substituting  the  far  field 

-i-Sv 


approximation  for  VoG(r|r“)  in  (2.10)  the  Kirchhoff  approximation  for 


the  scattered  acoustic  potential  is 


ip*(r)  -  -||^C(l+R)5f  -  (l-R)5i>i(?8)G(?|?8)-5dAc 


if--  -i  expCik(rg+r^)] 

"  Air  f  M*n  g^n ' - 5 — 5 — dA- 


i  f 
r  r 
s  s 


(2.18) 


in  which  r i  and  rf  are  distances  from  a  point  on  the  surface  A  and  the 

8  8  O 


source  and  observation  points  respectivelyt 
and 


M  -  Rv  -  p  ,  p-k(ni+  n*)  *  pa  +  p  a  +  pa 

xx  y  y  *  z 


v  ■  k(n^-n^)  •  v  p  +  v  a  +  va 
xx  y  y  z  z 


(2.19) 


The  results  of  the  analysis  for  the  non-depolariaed  scattered  electro¬ 
magnetic  field  is  the  same  as  (2.18)  with  the  exception  that  the  acoustic 
reflection  coefficient  is  replaced  by  the  Fresnel  reflection  coefficients 
for  horizonatally  and  vertically  polarized  waves  (Beckmann  and  Spizzi- 
chino,1963) . 
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2.4  The  Physical  Optics  Solution  and  the  So-Called  "Edge  Effect" 

In  order  to  apply  Green's  theorem  (2.9)  to  surfaces  that  are  not 
closed,  Beckmann  and  Spizzichlno  (1963)  suggest  that  "one  may  regard 
the  rough  surface  S  as  part  of  any  closed  surface  S'  of  which  only  S 
has  a  non-vanishing  reflection  coefficient,  the  remaining  integral 
over  S'-S  is  then  easily  shown  to  vanish  so  that  there  is  no  need  to 
Introduce  the  surface  S'". 

Furthermore  Beckmann  assumes  that  the  surface  S  is  such  that  the 

distances  r*  and  rf  from  the  source  and  the  observation  points  to  points 
s  s 

on  the  surface  S  are  much  larger  than  distances  between  two  points  on 
the  rough  surface.  In  this  case  it  is  convenient  to  define  a  new  coordinate 
system  rs(xs»ys»zs)  associated  with  the  surface  S  whose  origin  is  shifted 
from  the  original  coordinate  system  r(x,y,z)  such  that 


r(x,y,z)  -  r*+  zs(x8.ys»z8>  (2.20a) 

The  surface  S  is  defined  by 

f(xs*ys*2s)  “  zs  '  «<Vy8>  ■  °*  _X  <  XS  *  X*"Y  <  yB  <  Y  (2.20b) 

In  (2.20)  r^  is  the  constant  vector  joining  the  origins  of  the  two 

coordinate  systems  and  zg  ■  0  is  the  mean  plane  of  the  surface  S,  thus 

|C(*  »y_)d*  dy-  “  0  (2.21) 

J  8  S  S  8 

8  -i  -f 

In  this  case  for  points  on  the  surface  S,n  and  n  are  constant  and 

i  i  -i  -  f  f  -f  -  ,, 

r=r+n*r  r  =  r  -  n  *r  (2.22) 

so  s  so  s 

The  distances  from  the  source  and  observation  point  to  the  origin  of 

the  new  coordinate  system  are  r*  and  r^  respectively,  (see  Fig.  2.2). 

1  **s 

Thus  for  points  on  the  surface  S,  the  Incident  field  ^  (rQ)  is 


ikr  , -i  -  ,, -i  - 

.  .o  ikn  *r  .  ikn  *r 

,i,-Sv  ,-iv  e  s  „i  s 

(rQ)  -  g(n  )  - —  e  “  Eo  e 


(2.23) 


Figure  2.2  The  scattered  field  from  a  rough  surface. 


For  a  uniform  plane  wave  excitation  aaaumed  by  Beckmann  and  Spizzichino 
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(1963)  gCn1)  is  assumed  to  be  constant  and 


E1  -  1 
o 


(2.24) 


The  unit  vector  n  normal  to  S  is 

s  -  «/ivt  i  +  :.>/[(H>2  +  <!f>2  +  1]l1  <2-25) 

and 

dA  ■  dx  dy  /n*a  •  dx  dy  | Vf )  (2.26) 

O  8  S  Z  8  8 

Thus  subject  to  all  the  assumptions  made  by  Beclonann  the  scattered  field 
(2.18)  reduces  to  (Beckmann  and  Spizzichino ,1963;  Ishimaru,  1978) 


* a(x ) 


ikr 

i£~r  Is' 

4irr  1 


iv*r 


Vf  e 


s 


dx  dy 
s  Je 


(2.27) 


o  S 


where  the  integration  is  over  the  open  surface  S  which  has  a  non-vanishing 
reflection  coefficient.  For  a  flat  surface  normal  to  the  unit  vector 


n  *r 
o  s 


n«na+na+  n  _a 
o  ox  x  oy  y  oz  z 

n  *r  n  x  +  n  y 

^  os  _  .  ox  s  oy  s 

“  0  ,  f(zj  -  — - *  + - - - 1 — 


(2.28a) 

(2.28b) 


oz 


oz 


the  scattered  wave  is  given  by 
ikrf 


-  — ~  M-Vf4XY  sincC  (v%-  v,,  "^/n^JOsincC  (vy-Vj,  nQy/noz)Y] 


47Tr 


oy 


(2.29) 


in  which  sinca  ■  (sina)/a,  and 


M’Vf  -  (Rv  -  p)' 


(2.30) 


oz 


For  rough  surfaces  Beckmann  obtains  on  assuming  R  “  const,  and  on  integrating 
by  parts 


ikr 


r(r)  - 


ie 


47rr 


M*v'1 


iv*r  M 

e  dx8dys+  i 


riv-?  ?Y 

8 
e 


)  L 


dx 


M 

z 


' 

"  iv*r  " 

X 

_ i 

e  8 

dy 

. 

- 1 

*1 

(2.31) 
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Using  the  full  wave  approach  (Bahar ,1981a) ,  it  can  be  shown  that 
for  high  frequencies  R  in  (2.31)  must  be  evaluated  at  the  specular  point 
where  n  ■  v/v. 

The  last  two  terms  in  (2.31)  which  do  not  depend  on  the  shape  of  the 
surface  are  identified  by  Beckmann  as  an  "edge  effect."  These  terms  are 
considered  to  be  very  small  compared  to  the  first  term  in  (2.31)  and 
therefore  negligible  (Beckmann  and  Spizzichino,1963) .  Since 


2  2  2 

-  -  2  R(v  +  v  +  v  ) 

M*v  Rv  x  y  z 


(2.32) 


the  scattered  field  without  the  "edge  effect"  is  given  by 


ikr  ,  . 

.  o  _  2  i 
s  ie  Rv 

*  <r)  ‘  — T  —  • 

4ffr  z 

0  S 


dx  dy 
s  s 


(2.33) 


On  applying  (2.33)  to  flat  surfaces  (2.28)  one  gets 


o 


[“■  4XYsinc[(v  -v  n  /n  )X]sinc[(v  -v  n  /n  )Y]  (2.34) 
(  v  j  x  z  ox  oz  y  z  oy  oz 


The  results  for  scattering  by  flat  surfaces  of  area  A  (2.29) (which  con¬ 
tains  the  "edge  effect")  and  (2.34)  (which  does  not  contain  the  "edge 
effect")  are  now  compared  in  detail.  For  scattering  in  the  specular 
direction  with  respect  to  the  plane  ^0*^s  *  0  (2.28). 


-f_-i 

n  *  n  “  -n>  ■  -  v/v  ,  v  »  Ivl 

o  os  i-f  -ii  *  11 

ln  -n  I 

M*Vf  -  -(Rv-p) *v/vn  ■  -Rv^/vn  ■  Rv^/v 
oz  oz  z 


(2.35a) 


(2.35b) 


Thus  for  the  specular  case  the  "edge  effect"  vanishes  and  the  two  solu¬ 
tions  (2.29)  and  (2.34)  are  identical.  The  specularly  scattered  field 
from  the  flat  surface  n  *r  *0  (2.28)  is 

O  8 
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Cec(?)  "  ~~'7~  R  CO80s  A  (2‘36) 

ro 

In  which  R  is  the  reflection  coefficient  for  a  wave  incident  at  the 
angle  8g  upon  the  surface  nQ*r8  “  0«  Since  4  XY  is  the  projection  of 
the  surface  no*rg  *  0  on  the  xy  plane,  the  area  of  this  surface  is 


and 


A  - 


4  XY  4  XY 


n*a 


oz 


(2.37) 


v  -*■  2k  cos8  ,  v  -*■  vn  (2.38) 

a  z  z 

in  which  8g  is  the  angle  of  incidence  and  scatter  for  a  specularly 

oriented  surface,  n  m  n  ,  (2.35a).  Thus,  the  Physical  Optics  solu- 
o  os 

tions  (2.29)  and  (2.34)  for  specular  scattering  by  a  plane  of  arbitrary 
slope  satisfy  the  realizability  relationship  in  electromagnetic  theory. 

In  the  general  case  however,  they  could  differ  very  significantly.  To 
test  these  solutions  for  reciprocity  interchange  the  location  of  trans¬ 
mitter  and  receiver.  This  amounts  to  the  following  substitutions  in 


(2.27)  and  (2.33) 


h*  -n*  ,  n^+  -n*  ,  v  -*•  v  ,  p  ■*  -p  (2.39) 

Clearly  therefore  (2.33)  satisfies  reciprocity  while  (2.27)  does  not. 

Furthermore  on  setting  R  0  solution  (2.33)  vanishes,  as  must  be  the 
case  for  a  non-reflecting  surface.  However,  on  setting  R  -►  0  the  solu¬ 
tion  (2.27)  which  includes  the  so-called  "edge  effect"  reduces  to 
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dx  dy 
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(2.40) 


On  Integrating  over  the  surface  no*rg  ■  0  (2.40)  becomes 
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The  above  expression  vanishes  only  for  backseatter  (n£  “  -n*)  or  the 

specular  case  in  which  n  >  n  (2.35a)  and 

o  os 

-(n1*  nf)-nos  -  ^  -  0  (2.42) 

Thus  while  the  "edge  effect"  is  not  neglected  in  (2.27),  this  Physical 
Optics  solution  for  the  scattered  field  does  not  satisfy  reciprocity  nor 
does  it  vanish  for  the  non-specular  case  as  the  reflection  coefficient 
becomes  vanishingly  small.  On  the  other  hand  the  Physical  Optics  solution 
(2.33)  without  the  term  associated  with  the  "edge  effect"  satisfies  the 
reciprocity  relationship  (2.39)  and  vanishes  as  R  4  0.  It  also  satisfies 


the  realizability  relationship  (2.36).  However,  it  has  been  shown  that 
in  general  the  so-called  "edge  term"  is  not  necessarily  very  small  in 
comparison  to  the  first  term  in  (2.31)  which  corresponds  to  the  solution 
(2.33)  (Fung,  Moore  and  Parkins, 1965) .  In  these  cases  it  has  been  suggested 
that  "the  'edge  effect'  .  .  .  may  become  important  when  considering 
scattering  in  directions  far  off  specular".  Yet  it  has  been  shown  here 
that  on  retaining  the  "edge  effect",  the  Physical  Optics  solution  does 
not  satisfy  reciprocity  and  results  in  a  peculiar  non-vanishing  scattered 
field  for  the  non-specular  direction  (2.41). 

The  questions  that  need  to  be  resolved  therefore,  are  whether  the 
"edge  effect"  should  be  suppressed  only  when  it  can  be  demonstrated  that 
it  is  relatively  small  compared  to  the  scattered  field  and  if  it  needs 
to  be  included  in  the  solution  when  it  is  not  relatively  small,  how  does 
one  account  for  the  fact  that  the  resulting  solution  (2.27)  does  not 
satisfy  reciprocity  and  does  not  vanish  as  the  reflection  coefficient 
becomes  vanishingly  small? 
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2.5  Evaluation  of  the  Closed  Surface  Integral  for  the  Scattered  Field 
In  order  to  resolve  the  questions  raised  In  Section  2.4 ,  we  examine 
the  suggestion  that  "one  may  regard  the  rough  surface  S  as  part  of  any 
'closed  surface  S'  of  which  only  S  has  a  non-vanishing  reflection 
coefficient,  the  remaining  integral  over  S'-S  is  then  easily  shown  to 
vanish  so  that  there  is  no  need  to  introduce  the  surface  S"'.  To  this 
end,  the  starting  point  of  our  present  analysis  is  the  "closed  surface" 
integral  (2.18)  in  which  the  vectors  M,v  and  p  as  well  as  the  distances 
r*  and  r^  are  functions  of  position  on  the  surface  Aq  (spherical  waves). 

Divide  the  closed  surface  Aq  into  two  parts  and  and  let  the 
equation  for  the  surfaces  and  be  (see  Fig.  2.1) 

fx (v ,y  ,z)  -  z  -  C1(x,y)  -  0  (2.43) 


and 


f2(x,y,z)  -  z  -  £2(x,y)  -  0  (2.44) 

respectively,  such  that  on  S^,  n*az  >  0  and  on  S2»  n*az  <  0  and  on  the 

common  closed  line  boundary  C,  n*a  ■  0. 

z 

Thus 

K  >  0 


ndA  ■  dx  dy  J 
°  s  jb  • 

-Vf0  n*a  <  0  (2.45) 

l  2  z  " 

where  n  is  the  inward  normal  to  the  closed  surface  A  and 

o 

f1(x,y,z)  -  f2(x,y,z)  on  C  (n*az  -  0)  (2.46) 

Assume  that  the  reflection  coefficient  R  does  not  vary  rapidly  over  the 
surface  Aq  (i.e.,  ignore  scattering  due  to  changes  in  the  reflection 
coefficient)  and  assume  that  the  surface  Aq  is  not  excited  by  the  source 
at  grazing  angles  such  that 

|g(ni)vz|  <  K  ,  K  positive  real  const.  (2.47) 

Thus  retaining  only  the  far  field  terms  it  can  be  shown  that 
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4I  (Mg (n1) exp  ik(r*+  r*)/r*  r[  Vf)  =  V#.f 


■  i  F • (v  -v  Vf ) 

X 

(2.48) 

in  which 

rl  • 

(2.49a) 

aDd  rf  -  [(x  -  x)2  +  (y  -  y)2  +  (C  -  e)2]* 

XX  X 

(2.49b) 

and  f  and  £  are  defined  in  (2.43).  Denoting  the  value  of  f  (2.48)  on 
and  by  F^  and  respectively  and  using  the  divergence  theorem  in 
two  dimensions,  it  follows  that  since  on  C 

|v- (Fx-  F2)dx8dyg  -  Ij  (Fx-  F2)-(;xdy8-  iydxg)  -  0 

-  ij[F1-(v-vz7f)1  -  F2*(v  -vzVf)2]dxsdy8  (2.50) 

In  view  of  (2.45),  (2.50)  reduces  to 


j[V^lVVVVz2Vf2]dJtsdy8  5  f  vzf,SdAo  “  ^8(?> 

-  j[(F‘v)1  -  (5r*v)2]dxgdy8 

_  _  _  _  _  2 

and  since  M\v  -  (Rv  -p)*v  ■  Rv  ,  (2.18)  reduces  to 
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S  8 


dx  dy 

8  S 


(2.51) 


(2.52) 


dA_  -  dx  dy  /n  *a 
o  s  's  o  z 


Thus  the  above  expression  does  not  contain  the  term  associated  with  the 
"edge  effect."  The  integrand  In  (2.52)  corresponds  to  the  specularly 
reflected  wave  from  an  elementary  surface  of  area 

(2.53) 

Thus  the  integrand  of  (2.52)  is  the  contribution  to  the  scattered  field 
from  an  elementary  surface  and  may  be  regarded  as  proportional  to  the 
"specific  reflectance."  In  (2.52)  instead  of  integrating  over  the  closed 
surface  Aq  one  may  integrate  only  over  the  surface  S  for  which  R  i*  0. 
Furthermore,  if  the  distances  from  the  source  and  observation  points 
to  points  on  the  surface  are  very  large  compared  to  distances  between 
any  two  points  on  the  surface  S,  n1  and  n*  can  be  assumed  to  be  constant 
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(see  Section  2.4)  and  (2.52)  reduces  to 

s  -  iRv2  -i  Clk(ro+  ( 

♦  <*)  •  g(n  )exp - j — 2 - j  exp(ivrg)dxgdyg  (2.54) 

2  ro  ro  S 

Thus  the  terms  identified  with  the  "edge  effect"  (2.31)  should  be 
suppressed  even  when  they  are  not  relatively  small.  (Beckmann  and 
Splzzichino,1963) . 

It  should  be  noted  that  the  above  expression  was  derived  under 
the  assumption  that  |g(n*)/v2|  <  K  (2.47).  Thus  (2.53)  cannot  be  used 
for  excitation  at  grazing  angles.  For  grazing  angles,  the  Physical  Optics 
approach  is  not  valid  and  a  spectral  approach  must  be  used  (Bahar ,1981a) . 

To  account  for  shadowing  in  the  high  frequency  limit  one  should  note 

s  1.  if 

that  ip  -  ~\p  on  the  shadow  boundary  since  the  total  field  i|i  *  \p  +  \p 

vanishes.  Furthermore,  when  the  observation  point  is  in  the  shadow  region, 

g 

the  scattered  field  consists  of  the  "reflected  wave"  4>  and  the  "shadow¬ 
forming"  wave  (Morse  and  Feshbach,1953) . 

2 . 6  Cone lud ing  Remar  ks 

It  is  shown  here  that  the  discrepancies  between  the  different 
formulations  for  the  scattered  field  and  the  problems  associated  with 
the  so-called  "edge  effect"  in  the  Physical  Optics  solution  are  due  to 
the  premature  truncation  of  the  closed  surface  integral  expression  for 
the  scattered  fields.  On  applying  the  divergence  theorem  (in  two  dimensions) 
to  the  closed  surface  integral  (Section  2.5),  the  term  identified  as  the 
"edge  effect"  is  shown  to  vanish  identically  for  all  scatter  directions. 

Thus  using  the  Physical  Optics  approach  (which  cannot  account  for  wave 
diffraction  by  edges),  the  ao-called  "edge  effect"  term  must  be  sur- 
pressed  even  when  it  is  not  very  small  compared  to  the  scattered  field. 
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As  a  result  the  final  form  for  the  Physical  Optics  solution  (2.52) 
satisfies  reciprocity  and  the  realizability  relationships  and  it  can  be 
used  for  arbitrary  excitation.  The  integrand  in  (2.52)  is  identified 
with  the  "specific  reflectance"  (per  unit  area)  and  it  vanishes  as  the 
reflection  coefficient  for  the  rough  surface  vanishes. 
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2. A  Appendix 

For  the  convenience  of  the  reader,  the  formulations  for  the  scattered 
fields  derived  by  Eckart  (1953),  Parkins  (1967)  and  Boyd -Deavenpor t  (1973) 
are  summarized  here  for  the  Direchlet  condition  and  point  source  excita¬ 
tion  (s(r)  *  6(r)  (2.1)).  Thus  on  substituting 


ikr. 


<M?“)  -  / GaQ )  +  <J>8(?J)  -  0,  4»1(?*) 


,  (R  -1) 


(2A.1) 


and  the  Kirchhoff  approximation 


[H "  IS" + U- =  2  ^  ]-*  “ 2ik  *i,i: 


(2A.2) 


into  (2.10),  and  on  replacing  the  closed  surface  Aq  by  the  open  surface  S 
one  obtains  Parkins  formulation 


j  s1-:  ♦‘(fj)  »- 


ikr 


dS 


(2A.3) 


If  on  the  other  hand  (2A.1)  and  (2A.2)  are  substituted  into  (2.16a),  one 
obtains  Eckart* 8  formulation 


ikr; 

J  (ni-nf)*n  ^(r®)  j-  ds 
'  °  r 


(2A.4) 


Boyd  and  Deavenport  assume  that  the  Green's  function  G(r|rQ)  satisfies 
(2.7)  as  well  as  the  boundary  conditions  (2A.1)  and  (2A.2)  at  S.  Thus 
using  (2.16a)  as  the  starting  point  and  substituting  G(r|r®)  *  0, 
3G(r|r®)/3n  a  -2ik  n**n  G(r|r")  one  obtains  the  Boyd -Deavenport  formulation 


*■<?)  -  #  f 


ikr1 
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dS 


(2A.5) 


The  three  formulations  (2A.3),  (2A.4)  and  (2A.5)  will  obviously  give 
different  results  except  for  backscatter  n“*  ~n*.  Mote  also  that  in  all 


f  "  ™  -ir  ' 
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the  above  formulations  the  integrands  are  equal  at  the  specular  (stationary 


phase)  points  where  n  -  nQg  (2.35a)  and 


(2A.6) 


Thus  the  above  three  physical  optics  formulations  (2A.3),  (2A.4)  and 
(2A.5)  can  be  shown  to  be  in  agreement,  provided  that  the  respective 
integrals  are  evaluated  using  the  stationary  phase  method  (Bahar, 1981a) . 
However,  the  discrepancies  between  the  above  three  formulations  increases 
as  non-specular  scattering  becomes  more  important.  When  the  most  signifi¬ 
cant  contributions  to  the  scattered  fields  do  not  come  from  specular 
points  (as  in  the  case  of  backscatter  near  grazing  Incidence)  the  Physical 
Optics  or  Kirchhoff  approximations  are  not  valid  and  other  methods  must 
be  used  to  determine  the  scattered  fields  (Bahar,  1981a) . 

In  Section  2.5,  it  is  shown  that  either  (2.10)  or  (2.16a)  can  be  used 
for  the  starting  point  of  the  analysis  provided  that  Integration  is  over 

the  closed  surface  A  .  Thus  on  following  the  analytical  procedures  of 
■  o 

Section  2.4  and  noting  that 


.  ,  .  ,-i  -f,  -  2 

,  -i  -  ,  -f  -  k(n  -  n  )'v  v 

kn  -v  -  -kn  •  v  -  — * — - ‘ —  ■  — r  ■ 


,2n  -£  -i. 
k  (1-n  -n  ) 


(2A.7) 


all  the  above  three  formulations  can  be  shown  to  reduce  to  (2.5)  with 
R  ■  -1  and  g(n*")  “  1  prov’ded  that  in  (2A.3),  (2A.4)  and  (2A.5),  the 
integration  is  over  the  closed  surface  Aq. 
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3.0  DEPOLARIZATION  OF  THE  SCATTERED  RADIATION  FIELDS  BY  CONDUCTING 
OBJECTS  OF  IRREGULAR  SHAPE  ABOVE  ROUGH  LAND  AND  SEA-FULL  WAVE 
SOLUTIONS 
3 . 1  Background 

In  this  section  the  full  wave  approach  to  problems  of  scattering 
by  rough  surfaces  has  been  applied  to  the  problem  of  depolarization  of 
the  scattered  radiation  fields  by  objects  of  finite  conductivity  and 
irregular  shape.  In  the  analysis  complete  expansions  are  employed, 
exact  boundary  conditions  are  imposed  and  a  variable  coordinate  system 
that  conforms  with  the  local  features  of  the  irregular  surface  is  used. 

The  full  wave  solutions  are  expressed  in  forms  that  can  be  readily  com¬ 
pared  with  earlier  solutions  and  they  can  be  used  to  reconcile  the 
differences  and  bridge  the  wide  gap  between  these  solutions.  Thus,  the 
full  wave  solutions  for  the  backscatter  cross-section  are  shown  to 
reduce  to  the  Physical  Optics  solutions  when  the  high  frequency, 
stationary  phase  approximations  are  used.  Similarly,  for  slightly 
rough  surfaces  the  full  wave  expressions  reduce  to  the  perturbational 
solutions  for  Che  backscatter  cross-section. 

The  full  wave  solutions  are  shown  to  be  consistent  with  the  duality 
reciprocity  and  realizability  relationships  in  electromagnetic  theory. 

These  solutions  are  invariant  to  coordinate  transformations.  Since 
upward  and  downward  scattering  are  considered  in  the  analysis,  multiple 
scattering  and  shadowing  effects  can  be  taken  into  account  in  a  self- 
consistent  manner.  Thus,  the  total  scattered  field,  varies  continuously 
as  the  observer  moves  across  a  shadow  boundary  and  there  is  no  need  to 

i 

1 
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introduce  transition  terms  derived  from  other  theories.  The  full  wave 


approach  can  be  applied  to  deterministic,  periodic  rough  surfaces.  It 
can  also  be  used  to  determine  the  scattering  by  finite  scatterers  in 
the  presence  of  rough  land  or  sea. 

Scattering  by  random  rough  surfaces  is  dealt  with  in  detail  in  Sec¬ 
tion  3.5.  The  deterministic  scattering  problem  considered  here  in 
detail  provides  the  basis  for  such  an  analysis. 

3.2  Discussion 

In  this  paper,  the  full  wave  approach  has  been  applied  to  the 
problem  of  depolarization  of  the  scattered  radiation  fields  by  objects 
of  finite  conductivity  and  irregular  shape.  Since  the  full  wave  approach 
(Bahar ,1981a)  can  account  for  multiple  scattering,  the  irregularly  shaped 
objects  may  be  in  the  vicinity  of  rough  land  and  sea.  Thus  these  solutions 
could  also  be  used  to  distinguish  between  the  radar  returns  from  the 
irregularly  shaped  object  and  the  clutter  from  rough  surfaces  in  the 
background . 

The  principal  elements  of  the  full  wave  approach  are  as  follows: 

(Bahar ,1981a)  Complete  expansions  of  the  electromagnetic  fields  are 
employed.  Thus,  the  electric  and  magnetic  fields  are  expressed  in  terms 
of  the  radiation  term  as  well  as  the  lateral  wave  and  the  surface  wave 
terms.  Since  in  general  the  irregularly  shaped  objects  and  rough  surfaces 
depolarize  the  incident  waves,  the  complete  expansions  Include  both 
vertically  and  horizontally  polarized  waves  (Bahar, 1973a) .  Exact  boundary 
conditions  for  the  total  fields  are  imposed  at  the  irregularly  shaped 
surfaces.  Furthermore,  the  medium  is  characterized  by  the  electromagnetic 
parameters,  e  and  y  (which  may  vary  along  the  propagation  path)  and  it 
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is  not  necessary  Co  employ  approximate  impedance  boundary  conditions 
(Bahar  and  Rajan, 1979).  Precise  mathematical  procedures  are  used  in 
the  analysis.  Since  the  characteristic  functions  used  in  the  complete 
field  expansions  do  not  individually  satisfy  the  boundary  conditions  at 
the  irregular  surfaces,  in  general  the  field  expansions  do  not  converge 
uniformly  on  the  boundaries.  Thus,  for  example  in  order  to  avoid  term- 
by-term  differentiation  of  the  field  expansions,  use  is  made  of  Green's 
theorems  (Bahar, 1973b) .  Maxwell's  equations  for  the  transverse  compon¬ 
ents  of  the  electric  and  magnetic  fields  are  converted  into  a  rigorous 
set  of  generalized  telegraphists'  equations  for  the  forward  and  back¬ 
ward  traveling  wave  amplitudes.  These  first  order,  ordinary,  coupled 
differential  equations  can  be  solved  using  numerical  or  analytical 
techniques  (Bahar , 1973b) .  The  solutions  for  the  coupled  wave  amplitudes 
are  substituted  into  the  complete  field  expansions  and  the  steepest 
descent  method  is  used  to  obtain  the  expressions  for  the  scattered 
radiation  fields  (Bahar, 1981a) .  A  variable  coordinate  system  that 
conforms  with  the  local  features  of  the  irregular  boundary  is  used  to 
remove  the  restrictions  on  the  gradient  of  the  Irregular  boundary  (Bahar, 
1980).  (See  Appendix  3. A) 

The  principal  advantages  of  the  full  wave  approach  are  as  follows: 
The  full  wave  solutions  are  valid  for  all  incident  and  scatter  angles 
including  grazing  angles,  Brewster  angles,  scattering  in  the  specular 
direction  and  backscatter.  Since  upward  and  downward  scattering  are 
considered  in  the  full  wave  analysis,  multiple  scattering  and  shadowing 
effects  can  be  taken  into  account  in  a  self-consistent  manner.  Moreover, 
the  total  scattered  fields  vary  continuously  as  the  observer  moves  across 
a  shadow  boundary  and  there  is  no  need  to  introduce  transition  terms 
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derived  from  other  theories.  Since  exact  boundary  conditions  are  imposed 
in  the  analysis  and  there  are  no  restrictions  on  the  gradient  of  the  rough 
surface,  it  is  not  necessary  to  characterize  the  rough  surface  by  an 
approximate  Impedance  boundary.  Furthermore,  since  precise  mathematical 
procedures  are  followed  in  the  derivation  of  the  generalized  telegraphists’ 
equations,  the  full  wave  approach  can  be  used  even  when  the  characteristic 
functions,  used  in  the  complete  expansions,  do  not  individually  satisfy 
the  boundary  conditions  at  the  irregular  surfaces.  The  full  wave  solutions 
are  shown  to  be  consistent  with  duality,  reciprocity  and  realizability 
relations  in  electromagnetic  theory.  They  account  for  coupling  between 
the  radiation  fields,  the  lateral  waves  and  the  surface  wave  terms  of 
the  complete  field  expansions  (Bahar,  1977).  The  full  wave  solutions 
are  valid  from  low  frequencies  up  to  optical  frequencies  provided  that 
the  illuminated  surface  is  at  least  several  wavelengths  wide.  They  can 
be  expressed  in  a  form  that  can  be  readily  compared  with  earlier  solutions 
that  have  limited  applications.  Thus  these  solutions  can  be  used  by 
the  engineer  who  is  not  necessarily  familiar  with  the  analytical  tech¬ 
niques  used  in  the  derivations.  Furthermore,  the  full  wave  approach 
can  be  used  to  resolve  the  discrepancies  that  exist  between  earlier  solu¬ 
tions.  Thus  the  perturbations  solutions  (Rice,  1951;  Barrick,  1970), 
can  be  reconciled  with  the  corresponding  Physical  Optics  solutions 
(Beckmann  and  Splzzichlno,  1963;  Beckmann,  1968).  The  full  wave  approach 
can  be  applied  to  deterministic,  periodic  and  random  rough  surfaces. 

These  solutions  are  invariant  to  coordinate  transformations  and  they 
can  be  applied  to  problems  of  scattering  by  irregularly  shape*,  objects 
of  finite  conductivity  as  well  as  to  problems  of  scattering  by  rough 


land  and  sea. 
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In  Section  3.3,  the  problem  is  formulated  and  the  principal 
elements  of  the  full  wave  approach  are  summarized.  In  Section  3.4,  the 
full  wave  solutions  are  presented  for  the  scattered  radiation  fields  by 
irregular  objects  in  the  vicinity  of  rough  land  or  sea.  In  this  section 
the  effects  of  shadowing  and  multiple  scattering  are  accounted  for.  In 
Section  3.5  the  invariance  properties  of  the  full  wave  solutions  as 
well  as  the  duality  and  reciprocity  relationships  are  examined  in  detail. 
In  Section  3.6  the  stationary  phase  approach  is  used  to  obtain  the  high 
frequency  approximations  of  the  full  wave  solutions.  These  solutions 
are  compared  with  earlier  Physical  Optics  solutions  that  are  based  on 
the  Kirchhoff  approach.  In  Section  3.7  realizability  is  examined  and 
energy  conservation  is  shown  to  be  satisfied.  The  full  wave  expressions 
for  the  backscatter  cross-sections  are  derived  and  they  are  shown  to 
be  in  agreement  with  the  corresponding  perturbatlonal  solutions  if  the 
gradient  of  the  rough  surface  is  assumed  to  be  small.  Thus,  the 
discrepancies  between  the  perturbatlonal  and  the  Physical  Optics  solu¬ 
tions  are  resolved  and  the  wide  gap  that  exists  between  them  is  bridged 
by  the  full  wave  solutions.  In  Section  3.8  the  full  wave  solutions  are 
applied  to  grazing,  specular  and  Brewster  angles.  Thus,  if  at  grazing 
angles,  stationary  phase  points  do  not  exist  on  the  rough  surface,  the 
Physical  Optics  solutions  fail  even  at  very  high  frequencies.  Special 
forms  of  the  full  wave  solutions  are  given  for  good  conducting  boundaries 

3.3  Formulation  of  the  Problem 

To  determine  the  reflection  and  transmission  of  electromagnetic 
waves  at  the  interface  between  two  isotropic  media  characterized  by 
different  permittivities,  c,  and  permeabilities,  M,  it  is  necessary  to 
decompose  the  incident  (unperturbed)  electromagnetic  radiation  (far) 
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fields  into  vertically  and  horizontally  polarized  waves.  When  the 
Interface  is  a  plane  normal  to  the  unit  vector  ay  for  instance,  the 
plane  of  incidence  for  a  wave  traveling  in  the  direction  n*  is  defined 
as  the  plane  normal  to  the  unit  vector  a^,  where 


iHi  "  (.n-x  S^/ln1*  ay| 

Since  the  radiation  fields  are  normal  to  n*,  the  Incident  electric 

fields  can  be  expressed  in  terms  of  the  orthogonal  pairs  of  unit 

vectors  a  and  a  as  follows: 
ni  V  i 

-i  „Vi  -  .  Hi  - 

E  "  E  aVi  *  E  *Hi 


(3.1) 


where 


*Vi  ’  *Hix  n 


(3.2) 

(3.3) 


The  corresponding  magnetic  field  components  for  the  vertically  and 
horizontally  polarized  waves  are 
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in  which  ti  -  (y/e)^  is  the  intrinsic  impedance.  Thus 

u1  _  ,„Vi  -  ,  „Hi  -  .  -i 
H  -  (H  ayi+  H  aH1)  x  n 


(3.4) 


(3.5) 


When  the  boundary  between  the  two  media  is  irregular,  as  in  the  case 
of  a  rough  surface  f(x,y,z)  ■  0,  it  is  convenient  to  decompose  the 
radiation  fields  into  vertically  and  horizontally  polarized  components 
if  and  if  with  respect  to  the  "local  tangent  plane."  This  plane  is 
perpendicular  to  the  varying  unit  vector  n  which  is  the  outward  normal 
from  the  scattering  object  (see  Fig.  3.1).  Thus  in  this  case  the 
Incident  rsdiation  fields  are  expressed  as 


=i  _  Vi  -n  Hi  -n 
E  En  aVi  +  n  aHi 


(3.6) 
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where 


®Hi  "  “I 


(3.7) 


The  relationship  between  the  respective  components  of  the  vertically 
and  horizontally  polarized  waves  is 
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the  angle  between  the  planes  of  incidence  normal  to  au.  and  aa . .  Thus 

HI  HI 


C *  and  S*  are  expressed  in  terms  of  the  dot  and  scalar  triple  products. 


cj  “  c08*1  "  and  Slm  8ln^1  "  8Hi8Si  nl]  (3'10> 


Since  the  matrix  T  in  (3.9)  is  Hermetian 
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The  invariance  of  the  vector  E  in  (3.2)  and  (3.7)  requires  that 
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In  a  similar  manner  the  radiation  field  E  ,  scattered  in  the  direction 


-f 


n  can  be  expressed  in  terms  of  vertically  and  horizontally  polarized 


components  with  respect  to  the  reference  plane  normal  to  a^  or  with 


respect  to  the  plane  normal  to  n.  Thus 
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The  relation  nip  between  the  respective  components  of  the  vertically 
and  horizontally  polarized  waves  is 
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in  which  and  are  the  cosines  and  ainas  of  the  angles  between  the 
planes  of  scatter  normal  to  a^  and  ajj^.  The  expresslona  for  the  unit 
vectors  a^,  a^ ,  a”{  and  a^  aa  well  as  for  and  are  obtained 
from  the  above  expreaslona  for  a^,  a^,  a^,  a^,  C^j  and  on 
replacing  n*  with  n^. 

The  full  wave  expression  for  the  radiation  field  scattered  by  an 
elementary  area  dA  normal  to  the  unit  vector  n  can  be  expressed  in  matrix 
form  as  follows,  (see  Fig.  3.1,  Bahar,  1981a). 

dGf  -  GoC^nTfF(nfn,iiin)Ti  exp[iko<nf-  ni)-rg]dA-n  G* 

i  G  dC(nf ,ni)Gi  (3.15) 

o 

The  constant  Gq  la  given  by 

G  ■  -ik  exp[-ik  r^]/2irr*  (3.16) 

o  o  o 

jj 

in  which  k  »  w(y  e  )  is  the  free  space  wave  number,  y  and  e  are  the 
o  o  o  o  o 

free  space  permeability  and  permittivity  respectively,  and  rf  is  the 
distance  from  the  origin  to  the  observation  point.  The  assumed  exp(iwt) 
time  dependence  is  factored  throughout  this  work.  The  transformation 
matrices  T*  and  T*  are  defined  in  (3.10)  and  (3.14).  The  elementary 
area  dA  is  given  by 


dA  ■  n  dx  dz/(n *a^)  (3.17) 

and  the  matrix  G*  for  the  incident  fields  (3.4)  depends  upon  the 
excitation.  Thus  for  electric  and  magnetic  dipoles  oriented  in  the  ay 
direction  (Bahar  and  Raj  an,  1979) 


cPi  -  *pl  -  n. 
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ik  sin0*  exp[-ik  r*] 
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(3.18) 


in  which  r*  is  the  distance  from  the  dipoles  to  the  origin  and  J  (Amp- 
meters)  and  M  (Volt-meters)  are  the  electric  and  magnetic  dipole  moments. 
The  angle  0*  is  the  angle  of  incidence  in  medium  0  with  respect  to  the 
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reference  plane  normal  to  ay.  Thus,  (see  Fig.  3.1) 


cos0*  =  C*  ■  -rti‘a 
o  o  \ 


(3.19) 
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The  angle  6q  is  the  angle  of  incidence  with  respect  to  local  tangent 


plane  normal  to  n.  Thus  (see  Fig.  3.1  and  Appendix  3. A) 


„in  _in 
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(3.20) 


In  (3. IS),  r  is  the  position  vector  from  the  origin  to  the  elementary 
s 


area  dA  on  the  boundary  of  the  scatterer.  The  surface  of  the  scatterer 
is  given  by 


f (x.y.z) 


jy-h^x.z)  ,  n-ay  >  0 
|y-h2(x,z)  ,  n-ay  <  0 


(3.21a) 


Thus 


rs  “  X  \  +  h<x»z)®y  +  2  az  "  r  +  f(x,y,z)ay 


(3.21b) 

The  unit  vector  n  is  the  outward  normal  from  the  scatterer.  Thus 


n  -  ±  Vf/|Vf|  =  siny  cos6  a  +  cosy  a  +  siny  sinfi  a 

x  y  z 

-  ±  [-  h  a  +  a  -h  a  ]/(h2  +  1  +  h2)**  (3.22a) 

xx  y  z  z  x  z 

in  which  the  upper  and  lower  signs  are  used  for  surfaces  y-h^  or  y-h^ 
respectively  and 

h  -  9h/9x  and  h  “  9h/9z  (3.22b) 

x  z 

The  scatterer  is  characterized  by  the  complex  permittivity  and 
permeability  and  y^  respectively.  For  simplicity  in  this  paper  the 
scatterer  is  assumed  to  be  opaque,  therefore  transmission  through  the 
scattering  object  is  neglected.  The  unit  vectors  n*  and  n^  in  the 
directions  of  the  incident  and  scattered  waves  respectively,  can  be 
expressed  as  follows  in  terms  of  the  reference  coordinate  system 

n*  _  sinO*  c omd*  a  -  cos0*  a  +  sin0*  sind>*  a 
oTx  o  y  o  z 


(3.23) 
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and 


n^  ■  sin0^  cos$^  +  cos0^  a^  +  sin0^  ain$^  a^ 


(3.24) 


The  eleaenta  F1*^  of  the  2x2  acattering  matrix  FCn^.n111)  ,  (3.15)  are 

given  by  (Bahar,  1981a). 
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In  which  the  dimensionless  quantities  n^,  e^,  and  y^  are  n^  »  (e^y^/eQyo)  , 
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3.4  Full  Wave  Solutions  for  the  Scattered  Radiation  Field* — Effects  of 
Shadowing  and  Multiple  Scatter 

The  total  scattered  radiation  fields  are  obtained  by  Integrating 
the  expression  for  the  differential  scattered  field  (3.15)  over  the 
surface  of  the  scatterer  which  is  both  illuminated  and  visible  at  the 
observation  point.  All  four  elements  of  the  scattering  matrix  C*n  F, 
in  (3.25)  vanish  as 


Cdn  ■  -  n^’n  -  0 

0 

,  0*n  -*•  ♦  n/2 

(3.30a) 

or 

_fn  -f  -  A 

C  *■  n  *n  *  0 

,  6fn  t  tt/2 

0 

(3.30b) 

Thus,  the  scattered  radiation  fields  vanish  in  a  continuous  manner 
as  the  observer  moves  across  the  shadow  boundary.  In  the  shadow  regions 
however,  the  total  electromagnetic  fields  do  not  vanish  since  the  complete 
expansions  of  the  fields  constitute  the  radiation  term  as  well  as  the 
lateral  and  surface  wave  terms  (Bahar,  1973a, b,  1977,  1978). 

The  nonillumlnated  (shadow)  boundary  is  determined  by  points  on 
the  scatterer  rs  -  r^  that  satisfy  (see  Fig.  3.1) 

n^ntfr^)  -  g^x^.z^)  *  0  (3.31a) 

Thus  substituting  (3.22)  and  (3.23)  into  (3.31a) ,  the  boundary  of  the 

nonillumlnated  region  is  given  by 

tan8i[h*  cos^1  +  h*  sin^*]  *  -1  (3.31b) 

ox  z 

where 

hx  “  hx(fsl)  and  hz  “  V^l*  (3*31c) 

This  nonillumlnated  region  extends  to  points  on  the  scatterer,  rg  ■  r*2 

that  satisfy 
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where 


hl  “  h(?il}  *  h2  " 


(3.32c) 


Similarly  the  region  of  the  ecatterer  which  la  nonvieible  to  the  observer 

extends  from  points  on  the  acatterer  where  rg  •  r^  to  points  where 

-f 


-f 
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The  loci  of  rg  •  rg^  is  given  by 


5f-5(^i>  -  -  o 

(3.33a) 

Thus 

tan6^[h^  cos$^  +  h^  sin$*3  ■  1 
ox  z 

(3.33b) 

where 

>4  •  >  4-  >4^) 

(3.33c) 
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The  loci  of  r 

s 

■  r^  is  given  by 

<*s2  “  ^'"^sl*  “  ° 

(3.34a) 

Thus 

h2  -  ^  ‘  <*s2  -  Xsl)hx  +  {Zs2  *  4)hz 

(3.34b) 

where 

hj  -  h(^)  ,  h*  -  h(f^) 

(3.34c) 

Define  the  shadow  function  D(Fg)  such  that 


d(?8) 


(3.35) 


1  ,  illuminated  and  visible  region 
0  ,  nonilluminated  or  nonvisible  region 

where  the  nonilluminated  region  is  defined  by  (3.31)  and  (3.32)  while 
the  nonvisible  region  is  defined  by  (3.33)  and  (3.34).  Thus  using  (3.15) 
total  scattered  radiation  field  is  given  in  matrix  form  by 

Gf  -  G  [  ClnTfF  T1  expCik  (nf  -  n1) -r  ]D(r  )dA*n  G1 
oj  o  o  a  8 

A 

(3.36) 

in  which  the  matrix  T“F  T“  is  given  by 


=  Go  C(nf ,ni)Gi 
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(3.37) 
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For  elements  of  the  scattering  surface  where  ii,  ii*  and  ii*  are  coplanar 

([ii*  n*  n]  ■  0) ,  the  cross-polarlzatlon  terms  and  F1^  vanish.  If  in 

addition  n^,ni,n  and  ay  are  coplanar  and  «  S*  -  0,  the  matrix  T*  F  T* 

VH  HV 

becomes  diagonal  and  the  cross-polarization  terms  C  and  C  vanish 
(Bahar,  1980).  Note  however,  that  FP<fyo  and  CP<Vo  for  PfQ  even  if 
sin(<|>*-$*)  “  0. 

The  integrand  in  (3.36)  is  finite  as  long  as  the  area  of  the  scattering 
surface  (A  *  jdA*n)  is  finite.  Thus  the  expression  (3.36)  is  valid  even  as 
(n*a^)  «*  cosy  -*•  0  in  (3.17). 

Since  the  full  wave  analysis  accounts  for  upward  and  downward 
scattering  with  respect  to  the  reference  plane,  it  can  also  be  used  to 
determine  the  radiation  fields  due  to  multiple  scattering.  Thus,  if  the 
distance  between  regions  of  the  scatterer  (or  scatterers)  that  contribute 
to  multiple  scattering  is  sufficiently  large  to  justify  far  field 
techniques  (steepest  descent  or  stationary  phase  methods),  the  radiation 
field  initially  scattered  in  the  direction  n*  is  regarded  as  the  incident 
field  in  the  evaluation  of  the  multiply  scattered  field  (see  Fig.  3.2). 

Note  that  while  the  shadow  function  D  (r  )  ■  0  or  1  (3.35),  the 

8  S 

elements  of  the  scattering  matrix  F  vanish  in  a  continuous  manner  as 
the  observer  enters  the  shadow  region  (3.31)  and  the  total  radiation 
field  G*,  (3.36) , varies  continuously  as  observer  moves  across  the  shadow 
boundary. 

3.5  Invariance  Properties  of  the  Full  Wave  Solutions  and  the  Duality  and 
Reciprocity  Relationships 

In  order  to  show  that  the  full  wave  solution  for  the  scattered 
radiation  field  is  invariant  to  coordinate  transformation,  the  different 
coefficients  appearing  in  (3.36)  are  examined  in  groups.  The  Incremental 
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surface  area  of  the  scaCterer  is  dA  -  dA«n,  (3.17).  It  is  obviously 
invariant  to  coordinate  transformation. 

In  the  far  field  the  distance  from  the  source  to  points  on  the 
surface  of  the  scatterer  is 

Ir1^1  +  ?  |  -  r1  +  r  *5.  (3.38a) 

S  SI 

Thus  T*G*  exp[-ikQ  n*-^]  is  the  expression  for  the  incident 
(unperturbed)  far  field  at  the  surface  of  the  scatterer  decomposed  into 
vertically  and  horizontally  polarized  components  with  respect  to  the 
local  tangent  plane.  (See  Appendix  3. A)  It  is  invariant  to  coordinate 
transformation  in  view  of  (3.12). 

Similarly,  the  distance  from  points  on  the  surface  of  the  scatterer 
to  the  observation  point  in  the  far  field  is 

|r^  n^  -  r  )  -  r^  -  r  *n^  (3.38b) 

8  S 

Thus  the  expression  Gq  exp[lkQ  nf*rg]  is  also  invariant  to  coordinate 
transformation. 

Since  the  expressions  for  C^n  F ^  (3.25)  depend  upon  the  angles  of 
incidence  and  scatter  with  respect  to  the  local  tangent  plane  only, 
it  is  also  invariant  to  coordinate  transformation.  Finally,  in  view 
of  the  invariance  of  the  transformation  (3.14),  the  expression  for  the 
total  scattered  radiation  field  (3.36)  is  Invariant  under  coordinate 
transformation . 

Corresponding  to  the  transformations 

E  -*•  H,  H  ♦  -E,  J  M,  M  -3,  y  **  e  (3.39a) 

the  full  wave  solutions  (3.36)  can  be  shown  to  satisfy  the  following 
duality  relationships: 

FVV(e,v)  -  F^y.e)  ,  FHV(e,y)  -  -F^y.e)  (3.39b) 


r 


To  examine  the  reciprocity  relationships  the  locations  of  the  source 
and  observation  points  are  interchanged.  Thus,  under  the  transformation 


-i  _  -f  -1 
n  -n  ,  r 


(3.40) 


the  following  relationships  can  be  shown  to  be  satisfied 
T1^)  -  Tf(5f)  -  Tf(nf)T* 

Tf(-^)  -  TV)  -  ^(n1)1* 

F(-ni,-nf)  -  F(nf .n1)1* 


(3.41a) 

(3.41b) 

(3.41c) 


in  which  the  superscript  T  denotes  the  transpose  of  the  matrix  and  the 
symbol  *  denotes  that  the  nondiagonal  elements  of  the  matrix  are 
multiplied  by  -1.  Thus  on  noting  that  for  any  2x2  matrices  X  and  Y 

X*Y*  -  (XY)*  (3.42) 

the  following  relationship  for  C(n^,n*)*  (3.36)  is  obtained  (Bahar,  1981a), 


-i  -f.  r-.-f  -i,  -.T* 

C(-n  ,-n  )  •  [C(n  ,n  )] 


(3.43) 


and  the  total  scattered  radiation  fields  G  ,  (3.36)  satisfy  the  reciprocity 
relationships  in  electromagnetic  theory. 

3.6  High  Frequency  Approximations  of  the  Full  Have  Solution — 

Comparison  with  Physical  Optics  Solutions 

At  very  high  frequencies  the  major  contributions  to  the  scattered 
radiation  fields  (3.36)  come  only  from  regions  of  the  rough  surface  where 
the  phase  of  the  expression  exp[ikQrg* (n^-n*)]  is  stationary.  This  phase 
is  constant  when 


V[r  ‘(n^-n1)]  >V[(t  T  fa  )  •  (n^-n1)  ] 

8  y 

■  (n^-n*)  T  (n^-n^)*Sy  Vf  •  0 


(3.44) 


Thus  at  the  stationary  phase  points,  the  unit  vector  ng  (normal  to  the 
surface  f(x,y,z)  -  0  (3.22a),  is  in  the  direction  n*-n*  and 


now 
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-f  -1 
n  -n 

i-f  -it 
n  -n 


siny  cosd  a  +  cosy  a  +  siny  sinfi  a 

o  8  X  8  y  8  8  Z 


=  ±(-h  a  +  a  -  h  a  )cosy 
xs  x  y  zs  y  'a 


(3.45) 


The  incident  and  scatter  angles  0Q  with  respect  to  the  tangent  plane 
at  the  stationary  points  are 


.  4  ,  -f  -  r,  -f  -i r 

co.eo  -  -Co  •„,)  -  ■  L~t- J  ■ 


and 


l4CiCf-S1Sf 
o  O  0  0 

COS^^-**) 

2 

cos0^n8“  (n^*n  )  “  cosS^118  =  cos@S  =  Cs,-n*'n*  -  cos20S 
0  8  0  0  0  O 


(3. 


The  corresponding  angles  0*ns  and  0^ns  for  medium  are  obtained 


through  Snell's  law.  Thus 
,s  _  „s 


S’  -  S*/nr  and  C*  -  [l-tS*)2]** 


(3.47) 


Furthermore  at  the  stationary  points 


f  „  (n^x  n  ) * (n*x  n  ) 

,.fns  .ins.  s  s  , 

COS  (<|)  -<J)  )  -  - ; -  -  1 


.  ,.fns  .ins. 
sin(<()  -4>  ) 


|nfx  nj  If^x  ng| 

[ng  nf  n1] 

|nfx  SJ  l^x  ng| 


-  0 


cosy  ■  n  *a  ■  (C*  +  C^)/2  C8 
s  s  y  o  o  o 


sin6 


„f  .  .f  i  .  .1 
S  sia4>  -  S  sinA 
o  o 


,  CO86  " 


cf  .f  _i  . 1 
S  COS®  -  s  cos® 
o  T  o 


8  [(2  c8)2-  (C4+  C*)4]*  ’  '■''°U8  "  [(2  c8)2-  <?+  cf  ?]* 


(3.48a) 

(3.48b) 

(3.49a) 

(3.49b) 


1.8  f  8 

The  angles  ip  and  ip  are  given  by 

.  ,  Ci4Cf-Ci(l+cos  208)  .  .  Sf  sin(4»f-^i) 

_is  ,is  o  o  o  o'  _is  .  ...is  o 

Cj,  -  cos*  -  - j - ; -  ,  S.  -  sin*  -  - - - 

T  C1  .4.  9fi"  ▼  *  - 


S  sin  20 
o  o 


sin  20 


f-  f.  20!>  s„ 

_fs  .IS  O  O  O  O  _IS  .  .,18  O 

c  -  cosip  -  - S - - -  ,  -  sin*  -  - - - 

♦  Sf  sin  20s  ^  sin  20s 

o  o  o 

Thus  for  the  stationary  points  on  the  scatterer 


(3.50a) 


(3.50b) 


FPQ(nfn,nin)  -*■  0  PjlQ 


46a) 

46b) 


(3.51a) 


r 


CinFPP(nfn,nln)  C®R?  (C®)  2co829*RP  (0®)  ,  pg  _f  _± 

_2 -  —  -  C  kF®  F  (n  ,n  )  (3. 


Cf  +  C1 
o  o 


o  lo  o 


where  R?s  *  R?  (0®)  is  the  Fresnel  reflection  coefficient  for  vertically  or 
lo  10  o 

horizontally  polarized  waves  (P“V,H)  incident  at  an  angle  6®  upon  the 


plane  normal  to  n  .  The  function 
s 

2fts 

,  .  cos  9 

_  ,-f  -iv  o 


Fo(nf,ni)  - 


.  fe£+©1'1  fef-ei 

C1"  cos  ~  -r°  cos  — -j-  °- 
o  2  2 


l+CiCf-S1Sf  cos  (4>f  hJ)1) 
o  o  o  o  T 

c^c1-^) 

o  o  o 


(3.52) 


is  exactly  equal  to  the  function  F3 (0^ , 6£ , 4>f-«J>1>  derived  by  Beckmann  and 

Spizzichino  (1963)  using  a  Physical  Optics  approach  for  perfectly  conducting 

rough  surfaces.  Thus  if  there  are  points  on  the  illuminated  portion 

of  the  surface  of  the  scatterer  for  which  n  •  n  ,  the  Physical  Optics  limit 

8 

for  the  scattered  radiation  field  can  be  obtained  from  the  full  wave 
solutions  by  replacing  the  expression  for  C*n  Tf  F  T*ln  (3.36)  by  its  value 
at  the  stationary  points  and  by  replacing  the  shadow  function  D(r)  by 
unity.  Thus  in  the  Physical  Optics  limit  C^T*  F  T1  can  be  factored  out 
of  the  integral  (3.36)  and  the  expression  for  C  reduces  to 


C  -  CF  I 
o  o 


fs  Vs  is  _fs  Hs_is  fsVsisfsKsis 

RloS  “Sj>  *10^  lo^P  t  loSj» 


in  which  the  integral  1  is  given  by 

*(nf  ,n1,r-,A)  -  j  expCi^^-n1)*?^^*^ 

For  scatterers  that  are  very  good  conductors,  assuming  CQ  i*  n/2; 
R^*  -*■  1  and  R^*  -*■  -1,  (3.53a)  reduced  to 

cos(^f*-Hj>lB)  sindj^-Hj,**)' 

C  -  CiF  1  -  C*F  IU® 

o  o  ,  .  .  .  oo 

einO|r®-H|/18)  -coa^r*-^1  *) 


(3.53a) 


(3.53b) 


(3.54) 
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Thus  the  matrix  US  in  (3.54)  is  unitary. 

3.7  Realizability  and  the  Backacatter  Cross-Sections 

For  a  perfectly  conducting  plane  of  area  A  centered  at  the  origin 

and  normal  to  the  vector  n  ■  ii  ,  the  shadow  function  D(r  )  reduces  to 

s  s 

unity  and  (n^-n*)-rs  “  0.  Thus 

I(nf ,ni,ro,A)  -  A  n  *5  -  A  (3.55a) 

s  s  y  y 

in  which  A^  is  the  projection  of  the  area  A  on  the  reference  plane.  From 
(3.51b)  it  follows  that 


ch  - 


o  o  cosy 


-i  - 
-n  *n 


n  *a 
s  y 


(3.55b) 


For  this  case  the  scattering  matrix  C  (3.54)  is  therefore  given  by 

C(n^,n^)  “  A  cos08  U8 


(3.56) 


Since  IT  is  a  unitary  matrix  (3.54),  energy  is  conserved  independent  of 
che  angle  of  incidence  and  the  polarization  of  the  waves  incident  upon 
the  conducting  plane.  This  realizability  relationship  is  not  satisfied 
by  the  iterative  or  perturbat ional  solutions  (Bahar,  1981a). 

The  backscattered  radiation  fields  are  given  by  (3.36)  with  n^»  -n*. 


This  quantity  is  maximum  for  a  perfectly  conducting  plane  of  area  A 
centered  at  the  origin  and  normal  to  n  ■  n 


-n*  -  n^.  For  this  case  the 


backscattered  fields  are  obtained  from  (3.56)  on  substituting  0Q  -  0. 

f  a  ig 

Since  in  addition  <{/  +  ip  *0,  the  unitary  matrix  reduces  to  a 

diagonal  matrix.  Thus 


[cPQ(-s1,;1)]_„|-  *SpQ . 


PQ 

w 

where  is  the  Kronecker  delta.  Thus  the  maximum  backscattered  power 

2 

density  is  proportional  to  A  ,  the  square  of  the  area  of  the  perfectly 


(3.57) 
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conducting  plane  which  is  perpendicular  to  the  incident  wave  normal 
n*.  In  (3.57)  0q  is  the  angle  of  incidence  or  scatter  with  respect 
to  the  reference  plane  (y  ■  0). 

The  normalized  backscatter  cross  section  for  a  scatterer  of 
arbitrary  shape  is  defined  here  as 

0PQ-  [  |c^>Q(-ni,ni)  |/jD(rg)dA*n3^  (3.58) 

The  quantity  in  the  denominator  of  (3.58)  is  the  area  of  the  illuminated 

surface  of  the  scatterer.  Thus  the  maximum  value  of  the  normalized, 

PO 

dimensionless  quantity  o  defined  in  (3.58)  is  unity. 

When  the  Illuminated  portion  of  the  surface  of  the  scatterer  is 

slightly  rough  such  that  n  =  a^  the  perturbational  approach  is  valid. 

Thus  the  full  wave  solutions  for  the  backscattered  cross  sections  reduce 

to  the  perturbational  solutions  on  replacing  n  by  a^  in  the  expressions 

(3.25)  and  on  ignoring  the  shadow  function  D(r  )  (Barrick,  1970).  In 

& 

this  case  the  expression  for  C  (3.36)  simplifies  considerably  since 
C*nT^  F  Ti  beconies  independent  of  the  variable  of  Integration  and 

C  -  (cV  F  Ti)-^-^*I(nf ,ni,rg,A) 

=  c  KS^S1.?  ,A)  (3.59) 

\J  s 

where  the  integral  I  is  defined  by  (3.53b)  and  the  matrix  Cq  is  defined 
by  (3.59).  Similarly,  if  there  are  points  on  the  illuminated  portion 
of  the  surface  of  the  scatterer  for  which  n  ■  ng,  at  high  frequencies 
one  can  use  the  Physical  Optics  approximations  derived  in  Section  3.6. 

In  this  case  n  is  replaced  by  n  in  (3.36)  and  assuming  that  the  specular 
points  are  in  the  visible  portion  of  the  surface  of  the  scatterer,  D(rg) 
is  set  equal  to  unity.  Thus  for  the  high  frequency  limit  the  expression 
for  C  (3.36)  also  simplifies  significantly.  The  expression  for 
cJV  FT1  becomes  Independent  of  the  variable  of  integration  and 
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C  -  (C*  Tf  F  T1)—  »ni,rg,A) /n*a 

s  7 

=  I(nf,n1,rs,A)  (3.60) 

where  C  is  defined  by  (3.60). 

Thus  for  these  special  cases  (small  slope  or  high  frequency)  the  expression 
for  the  normalized  backscatter  cross  section  (3.58)  reduced  to 


tPQ 


in  which 


■  |  exp[-2ini*rs]dxdz |^/[|D(rg)dA*n]^ 

A  A 

|CP^|2  (  exp[-2ini,(r  -  r ')]dx'dz 'dxdz/[ |D(r  )dA*n]^ 
9>  iJ  S  8  J  8 


r'  "  x'a  +  h(x',y')a  +  z'a 

s  x  y  z 


(3.61) 


(3.62) 


PQ 

Since  the  full  wave  solutions  for  the  backscatter  cross  section  0  x  (3.58) 

PO 

reduce  to  the  perturbatlonal  form  Oq  for  slightly  rough  surfaces 

PO 

(Barrick,  1970)  and  to  the  Physical  Optics  form  om  when  the  stationary 
phase  approximations  are  used,  the  full  wave  solutions  provide  a  direct 
connection  between  these  special  forms  for  the  backscatter  cross  section 
and  reconcile  the  significant  differences  between  them.  Moreover,  the 
full  wave  method  can  be  applied  to  a  wide  class  of  problems  that  cannot 
be  solved  by  either  the  perturbatlonal  or  the  Physical  Optics  approach. 


3.8  Properties  of  the  Full  Wave  Solutions  at  Grazing,  Specular  and 

Brewster  Angles  and  Applications  to  Random  and  Periodic  Rough  Surfaces 

For  grazing  incident  or  scatter  angles,  with  respect  to  the  local 

tangent  planes,  (normal  to  the  vector  n)  all  the  full  wave  expressions 

for  C  nFP^(n^,ni)  (3.25)  vanish.  The  corresponding  expressions  in  the 
o 

perturbatlonal  solutions  vanish  only  for  grazing  incident  angles  with 
respect  to  the  reference  plane  (normal  to  a^).  On  the  other  hand,  the 
Physical  Optics  expression  for  (3.25)  vaniehs  only  for  grazing  angles 
with  respect  to  the  stationary  phase  plane  where  n  ■  ng.  The  property 
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of  full  wave  solutions  at  grazing  angles  with  respect  to  the  local 

tangent  plane  is  responsible  for  the  fact  that  the  full  wave  solutions 

for  the  scattered  radiation  fields  vanish  in  a  continuous  manner  as 

the  observer  moves  across  a  shadow  boundary.  For  points  on  the  Irregular 

boundary  where  the  local  tangent  plane  coincides  with  the  stationary 

phase  plane  normal  to  n  ■  ng  (locally  specular  scatter) ,  the  full  wave 
PO 

expressions  for  F  x  reduce  to 

FPP (nf n ,ivin)  -  RPP(e!)  ,  P-V.H} 

P*Q 


JPQ.-fn  -in.  A 
F  (n  ,n  )  •  0  , 


n  "  n 


(3.63) 


PP  PP  s 

The  corresponding  perturbational  solutions  for  F  reduce  to  R  (6  ) 

o 

only  for  specular  angles  with  respect  to  the  reference  plane  normal  to 

-  po 

a^.  Similarly,  the  perturbational  solutions  for  F  ^(P^Q)  vanish  only 

for  specular  angles  with  respect  to  the  fixed  reference  plane. 

For  the  specular  case  with  respect  to  the  reference  plane  the 

pp  pp  e 

Physical  Optics  solution  for  F  reduces  to  R  (0Q)  for  the  entire  rough 

surface  Including  those  portions  that  are  not  specularly  oriented  with 

respect  to  transmitter  and  receiver.  Furthermore,  while  the  Physical  Optics 

PO 

approach  does  not  provide  any  expression  corresponding  to  F  (P/Q) ,  the 

pn 

full  wave  expression  for  F  vanishes  only  at  the  specular  points  on  the 

irregular  boundary.  It  is  interesting  to  note  that  if  (the  incident 

or  scatter  angle  with  respect  to  the  stationary  phase  plane)  is  equal 

«B 

to  the  Brewster  angle  oq ,  the  Physical  Optics  solution  for  the 
scattered  radiation  field  vanishes.  Similarly,  if  the  Incident  on 
scatter  angles  with  respect  to  the  reference  plane  (normal  to  a^)  equals 
the  Brewster  angle,  the  perturbational  solution  also  vanishes.  The 
full  wave  solution  on  the  other  hand  does  not  vanish.  The  integrand 
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in  Che  expression  for  Che  full  wave  soluCions  (3.36)  vanishes  only  ac 

i  f  B 

Che  sCaCionary  phase  poinCs  where  6  >6  -  0  .  Thus.  Che  full  wave 

o  o  o 

soluCions  for  scacCering  by  rough  surfaces  does  noc  vanish  even  if 
Che  incidenC  and  scaCCer  angles  wich  respecC  Co  Che  stationary  phase 

plane  (normal  Co  n  )  or  wich  respecC  Co  Che  reference  plane  (normal 

8 

Co  a^)  are  equal  Co  Che  Brewscer  angle. 

For  very  good  conduce ing  boundaries  wich  Je^J  »  1  and 
yr  »  1  (|nrl  <<:  1)*  the  expressions  for  C^nF^  (3.25)  simplify  considerably. 

cJV™  -  2[cos(4>fn-<|>in)  -  SlnSfn]/(Cin+Cfn)  (3.64a) 

o  o  o  o  o 

-  -2C*nC*n  cos (4>f n-<t>in) / (C*n+C*n)  (3.64b) 

-  -2sin(^fn-4»itl)C^n/ (C^n+C^n)  (3.64c) 

C*nFVH  -  -2sin($f n-^in)C^n/ (C^n4C^n)  (3.64d) 

HH 

Wich  Che  excepCion  of  F  Che  above  expressions  cannoC  be  used  for  grazing 

angles  if  Che  boundary  is  highly  conducCing.  In  general,  the  approximate 

formulas  for  C^F^  (3.64)  must  be  restricted  not  only  to  good  conducCing 

boundaries  but  also  to  angles  of  incidence  or  scatter  that  are  less 

£ 

than  the  pseudo-Brewster  angles  6q.  For  vertically  polarized  waves 

Co  -  C®  -  COS0®  -  nrC®  -  [(n^-D/e^l]5*  -  (er+l)",s  +  e^.  (3.65) 

Thus  (3.64)  can  be  used  only  if  C*n  and  C^n  are  greater  than  (e^)"5*. 

For  perfectly  conducCing  boudaries,  the  steepest  descent  method, 
used  in  evaluating  Che  radiation  fields  at  grazing  angles,  needs  to 
account  for  the  poles  in  the  vicinity  of  the  saddle  points.  Thus  the 
apparent  singularity  in  the  expression  for  the  fields  at  grazing  angles 
over  perfectly  conducting  surfaces  is  removed  (Bahar,  1981b). 

The  full  wave  solutions  (3.36)  can  be  applied  to  periodic  structures 
by  multiplying  the  expressions  for  the  scattered  radiation  fields  due 


to  a  single  element  of  the  periodic  rough  surface  by  tvo  dimensional 
array  factors  (Bahar,  1980).  Implicit  in  these  array  factors  are  the 
expressions  for  the  grating  angles. 

The  expressions  for  the  full  wave  solutions  depend  upon  the  profile 

h(x,z)  of  the  rough  surface  and  upon  its  gradient.  Thus  in  order  to 

determine  the  statistical  average  of  (3.36)  and  its  variance, 

<CPQ>  and  D{CP^}  respectively,  it  is  necessary  in  general  to  know  the 

statistics  of  the  random  variables  h,  h  and  h  .  For  two  special 

X  z 

cases  (the  small  slope — perturbational  solution  and  the  stationary 
phase — geometrical  optics  solution)  considered  in  this  work,  the 
statistical  average  and  the  variance  of  the  scattered  radiation  fields 
can  be  expressed  in  terms  of  the  one  and  two  dimensional  characteristic 
functions  respectively  (Bahar,  1981a). 

3.9  Concluding  Remarks 

The  full  wave  approach  is  applied  in  this  paper  to  the  problem  of 
depolarization  of  the  scattered  radiation  field  by  an  object  of 
irregular  shape  and  finite  conductivity.  The  principal  elements  of 
the  full  wave  approach  are  outlined  in  the  introduction  (Section  3.3). 
The  full  wave  solutions  are  presented  in  a  form  that  can  be  readily 
compared  with  earlier  solutions  (Section  3. AX  They  are  used  to  resolve 
the  discrepancies  between  the  earlier  solutions  and  to  bridge  the  wide 
gap  that  exists  between  them.  Realizability,  reciprocity  and  duality 
relationships  in  electromagnetic  theory  are  examined  and  the  full  wave 
solutions  are  shown  to  be  Invariant  to  coordinate  transformations. 

The  full  wave  approach  can  also  be  used  to  determine  the  scattered 
radiation  fields  due  to  lateral  variations  in  the  permittivity  c  and 


permeability  y  (Bahar,  1973a, b).  Thus  they  can  be  applied  to  mixed  path 
propagation  problems. 

Since  the  full  wave  approach  has  been  generalized  to  inhomogeneous 
multilayered  structures  of  arbitrarily  varying  thickness,  it  can  also 
be  applied  to  problems  in  which  Irregularly  shaped  objects  are  imbedded 
in  the  earth’s  crust.  If  the  medium  of  the  scatterer  is  nondisslpatlve, 
transmission  through  the  scatterer  can  also  be  accounted  for  in  the 
analysis. 

When  the  transmitter  or  receiver  are  near  the  scatterer,  the  lateral 
wave  or  surface  wave  terms  of  the  full  wave  expansions  could  be 
significant  especially  near  a  shadow  region.  In  these  cases  the  full 
wave  approach  can  be  used  to  determine  the  coupling  between  the  radiation 
term  and  the  lateral  and  surface  wave  terms  (Bahar,  1977). 
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The  relationship  between  the  reference  coordinate  system  X(x,y,z)  and 
the  variable  local  coordinate  system  X(x,y,z)  associated  with  the  local 
tangent  plane  normal  to  n  (3.22)  can  be  expressed  as  follows: 


X  - 


Thus,  x  •  (sum  on  j  ■  1,2,3) 


r  i> 

I'l 

x 

X 

-2 

2 

X 

-  A 

X 

-3 
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lx  J 

lx 

(3A.1) 
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where  x  ■  x,  x  “y,  x  -  z,  etc.  The  elements  a^  of  the  3x3  trans¬ 
formation  matrix  A  are 


‘ij 


vaj 


i,j  “  1,2,3. 


(3A.2) 


where  a,  ■  a  ,  a„  *  a  and  a_  «  a  are  the  unit  vectors  of  the  reference 
1  x  2  y  3  z 

coordinate  system  and  n^  are  unit  vectors  normal  to  the  coordinate  surfaces 


of  the  local  coordinate  system.  Thus 

n1  ■  n  x  (a^x  n)/|a^x  n| 


n2  -  n 

i»3  -  (SjX  n)/|a1x  n| 


(3A.3) 

(3A.4) 

(3A.5) 


Thus  the  unit  vectors  n^,  ,  and  n^  are  orthogonal .  The  vector  n2  is 
normal  to  the  local  tangent  plane  (3.22)  while  n^  and  n^  lie  in  the  local 
tangent  plane.  For  a  horizontal  surface  h(x,z)  -  const.,  the  unit  vectors 
n^  are  equal  to  a^  and  A  reduces  to  the  identity  matrix.  In  general  the 
transformation  matrix  A  is  hermetian,  thus  the  determinant  of  A  is  equal 
to  unity  and 


(3A.6) 
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4.0  SCATTERING  CROSS  SECTIONS  FOR  COMPOSITE  RANDOM  ROUGH  SURFACES 

—FULL  WAVE  ANALYSIS 

4.1  Background 

In  this  work  the  full  wave  approach  to  rough  surface  scattering  is 
applied  to  composite  models  of  rough  surfaces.  It  is  shown  that  both 
specular  point  scattering  as  well  as  Bragg  scattering  are  accounted  for 
in  the  analysis  in  a  self-consistent  manner.  The  results  are  compared 
with  earlier  solutions  based  on  a  combination  of  Physical  Optics  and 
perturbation  theories.  Using  the  full  wave  approach  it  is  not  essential 
to  decompose  the  rough  surface  into  individual  surfaces  with  different 
roughness  scales  unless  it  is  desired  to  separate  the  specular  point 
contribution  from  the  Bragg  contribution  to  the  scattering  cross  sections. 
The  decomposition  of  the  rough  surface  not  only  enhances  one's  physical 
insight  but  also  facilitates  the  numerical  evaluation  of  the  scattering 
cross  sections.  Shadowing  is  also  accounted  for  in  the  analysis. 

4.2  Discussion 

Physical  Optics  and  perturbation  theories  have  been  applied  to 
problems  of  scattering  of  electromagnetic  waves  from  rough  surfaces 
f(x,y,z)*y  -h(x,z)*0  (Valenzuela,  1968) .  However,  these  theories  can 
only  be  applied  to  a  limited  class  of  rough  surfaces.  Thus  perturbation 
theory  (Rice,  1951;  Barrick,  1970)  can  be  applied  to  problems  in  which 
it  is  usually  assumed  that 

k2<h2>  «  1,  3h/3x  -  hx  «  1,  3h/3z  -  hz  «  1  (4.1) 

2 

in  which  <h  >  is  the  mean  square  of  the  rough  surface  height  and  kQ  is 
the  wave  number  of  the  electromagnetic  wave.  Physical  Optics,  which  is 
based  on  the  Kirchhoff  approximations  of  the  surface  fields  (Beckmann, 
1968),  is  applicable  to  surface  for  which  the  radii  of  curvature  of  the 
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rough  surface  are  large  compared  to  the  electromagnetic  wavelength  X. 

Wright  (1966)  and  Semyonov  (1966)  apply  the  above  theories  to 

composite  surfaces  made  up  of  irregularities  that  are  both  small  as  well 

as  large  compared  to  the  wavelength  X.  More  recently  Brown  (1978,  1980) 

applied  Physical  Optics  as  well  as  Burrow's  perturbation  theory  (1967) 

to  derive  the  expression  for  the  backscatter  cross  section  from  perfectly 

conducting  rough  surfaces  in  terms  of  a  sum  of  two  backseat terer  cross 

sections.  The  first  is  the  backscatter  cross  section  for  the  surface 

with  the  large  scale  roughness  h^,  while  the  second  is  the  backscatter 

cross  section  associated  with  the  small  scale  roughness  h^.  In  his  work. 

Brown  (1978,  1980)  assumes  that  the  radii  of  curvature  of  the  surface,  h^, 

2  2 

is  larger  than  the  wavelength  X  and  k  <h  >  >>  1.  In  addition,  he  assumes 

O  Aj 

that  hg  satisfies  the  conditions  (4.1).  Thus  in  his  work  the  specifica¬ 
tion  of  the  wavenumber  kj  (where  spectral  splitting  is  assumed  to  occur), 

2  2 

is  based  upon  the  characteristics  of  the  small  scale  structure  (k  <h  >  «  1) 

os 

rather  than  upon  the  characteristics  of  the  large  scale  surface  (Brown, 

1978;  Hagfors,  1966;  Tyler,  1978).  In  Brown's  analysis  the  backscatter 
cross  section  associated  with  the  surface  with  the  small  scale  rough¬ 
ness  is  expressed  in  terms  of  a  two  dimensional  convolution  of  transforms. 

In  the  composite  models  of  Wright  (1968),  Semyonov  (1966),  and  Valen¬ 
zuela  (1968)  which  are  "mostly  based  on  physical  considerations"  the 
rough  surface  is  approximated  by  "patches"  of  slightly  rough  surfaces 
that  ride  the  large  waves.  Thus  in  their  work  the  scattering  cross 
section  associated  with  the  surface  with  the  small  scale  roughness  is 
obtained  by  averaging  over  the  distribution  of  slopes  of  the  large 
scale  surface  roughness. 

In  this  work  the  full  wave  approach  to  rough  surface  scattering 
is  extended  to  composite  models  of  the  rough  surface  with  moderate 
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mean  square  slopes  (Bahar,  1981fe,c).  Since  the  full  wave  approach 
accounts  for  both  specular  point  scattering  as  well  as  Bragg  scattering 
in  a  self-consistent  manner  (Bahar,  1981b ,c)  without  introducing 
perturbation  and  Physical  Optics  theories,  it  is  not  necessary  to 
decompose  the  rough  surface  into  two  surfaces  with  large  and  small 
roughness  scales,  h^  and  hg  respectively.  Thus  the  need  to  specify  kj, 
the  wavenumber  where  spectral  splitting  is  assumed  to  occur,  does  not 
necessarily  arise  when  the  full  wave  approach  is  used.  However,  the 
decomposition  of  the  rough  surface  enhances  one's  physical  insight 
and  also  facilitates  the  numerical  evaluation  of  the  scattering  cross 
sections. 

In  Section  4.3  the  problem  is  formulated  in  terms  of  the  full 

wave  approach  and  the  principal  results  from  earlier  analysis  are 

summarized.  In  order  to  compare  the  full  wave  solutions  for  the 

scattering  cross  sections  with  earlier  solutions  (Valenzuela,  1968; 

Brown,  1978,  1980)  it  is  also  assumed  in  Section  4.4  that  the 

composite  surface  can  be  decomposed  into  two  statistically  independent 

surfaces  with  large  and  small  roughness  scales  h^  and  hg.  It  is  shown 

that  the  total  scattering  cross  section  can  be  written  as  a  weighted 

sum  of  individual  cross  sections.  The  cross  section  associated  with 

i  It  1 2 

specular  point  scattering  is  multiplied  by  | X  I  the  magni  ude  squared 
of  the  characteristic  function  for  the  small  scale  surface  height  h 

K 

(4.30a).  The  cross  section  associated  with  Bragg  scattering  is  expressed 
in  terms  of  an  Integral  over  the  slopes  of  the  rough  surface  (4.44)  as 
predicted  on  the  basis  of  physical  considerations  (Valenzuela,  1968). 
This  term  is  shown  to  be  in  agreement  with  perturbation  theory  (Rice, 
1951;  Barrick,  1970).  The  integrand  in  this  term  is  proportional  to 
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the  surface  height  spectral  density  VKv^.v^)  whose  arguments  are  the 
two  components  v_  and  v_  of  the  fector  v  ■  k  (n^-n*)  in  the  local  tangent 

X  2  ® 

plane  (n^  and  n*  are  unit  vectors  in  the  direction  of  the  scattered  and 
incident  wave  normals).  The  full  wave  solution  can  also  be  expressed 
in  terms  of  the  rough  surface  tilts  "in"  and  "perpendicular"  to  the 
plane  of  Incidence  and  thereby  compared  with  the  results  of  Valenzuela, 
et  al.  (Valenzuela,  Lang  and  Daley,  1971).  It  is  also  shown  that  since 
Brown  (1978,  1980)  assumes  that  the  small  scale  surface  height  correla¬ 
tion  function  is  a  function  of  distance  in  the  reference  (mean)  plane 
and  not  distances  along  the  large  scale  surface  as  in  this  work  and 
implicitly  in  Valenzuela’s  work,  the  full  wave  results  are  in  agreement 
with  Brown's  results  only  if  the  large  scale  surface  roughness  has  small 
mean  square  slopes.  Shadowing  is  accounted  for  explicitly  in  this  work 
as  in  the  work  by  Brown  (1978,  1980). 

Throughout  this  work  an  exp(iu)t)  time  dependence  is  assumed. 

4.3  Formulation  of  the  Problem 

For  the  incoherent,  diffuse  field,  the  scattering  cross  section  per 
unit  area  is  defined  as  (Ishlmaru,  1978), 

<0PQ>  -  47T(rf)2<|Epf-  <Epf>|2>/Ay|EQ1|2  ,  P,Q-V,H  (4.2) 

in  which  the  first  superscript  P«V,H  indicates  the  polarization  (vertical 
or  horizontal)  of  the  scattered  field,  EP^  and  the  second  superscript 
Q*V,H  indicates  the  polarization  of  the  incident  field,  E^.  The 
projection  of  the  area  of  the  rough  surface  A  on  the  reference  plane 
normal  to  a  is  A  and  rf  is  the  distance  from  the  origin  (associated 

y  y 

with  the  rough  surface)  and  the  observation  point  in  the  far  field  (see 
Fig.  4.1).  Using  the  full  wave  solutions  for  the  scattered  radiation 
fields  (Bahar,  1981a),  the  expression  for  the  scattering  cross  section 


Incident 

Planer 


Figure  4.1  Plane  of  Incidence,  scatter  plane  and  mean  (reference)  plane 
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for  random  rough  surfaces  is  given  by  (Bahar ,  1981b) 


<a 


pQ> 


k  r 

iH 


SPQ  exp[ivy(h-h')]- 


DP% 


<  ^  exp(lVyy)> 


n*a 


■  exp[iv  (x-x’)+iv  (z-z ')]dxdz  dx'dz *>  (4.3) 

x  z 

in  which  k  Is  the  free  space  wavenumber.  The  vector  v  in  the  cartesian 
o 

coordinate  system  is 


v  -  k  (n^-n*)  *  v  a  +  v  a  +  v  a 
o  xx  yy  zz 


(4.4) 


in  which  n*  and  n^  are  unit  vectors  in  the  direction  of  the  incident 


and  scattered  waves  respectively. 

The  characteristic  function  x(vy)  given  by 


X(vy> 


exp(ivyh)p (h)dh 


(4.5) 


in  which  p (h)  is  the  surface  height  probability  density  function. 
The  coefficient  SP^  is 


gPQ  _  PPQ(r)U(r)  DP^*(r’)U(r ') 


(n-ay) 


(n’-ay) 


p 


n*a 


(4.6) 


in  which  PP^(r)  depends  upon  the  polarization  (P,Q«V,H)  and  directions 
(n*  and  n^)  of  the  incident  and  scattered  waves,  the  unit  vector  n 
normal  to  the  rough  surface  and  the  electromagnetic  parameters  of  the 
medium  of  propagation  and  the  scatterer  (Bahar,  1981b).  The  shadow 
function  U(r)  is  unity  when  the  surface  is  both  illuminated  by  the  source 
and  visible  at  the  observation  point  and  zero  otherwise  (Sancer,  1969). 
The  symbol  <  >  denotes  the  statistical  average.  Thus  assuming  that  either 
the  rough  surface  radii  of  curvature  are  large  compared  to  a  wavelength 
or  the  rough  surface  slope  is  small  and  that  the  rough  surface  heights 
and  slopes  are  statistically  independent  (Bahar,  1981c), 
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<SPQ[exp  iVy(h-h')  -  !x(vy)|2>  -*• 


tZ\n tZ\  I  ^ 


f  DfV(r)U(r 
'  n-S 


[exptiv  (h-h ' )  ) - 1 X (Vv)  |  3p (S»0)W(h ,h *  )dhdh  'dh  dh  dU 


(A.7) 


The  joint  rough  surface  height  probability  density  function  is  W(h,h') 
and  the  joint  probability  density  function  p(n,U)  can  be  expressed  in  terms 
of  the  conditional  density  p(U|n)  (Sancer,  1969) 

p(n,U)  -  p(n)p(U|n)  (A. 8) 

in  which  p(n)  is  the  density  function  of  the  gradient  of  the  rough  surface 


V£  -  V(y-h(x,z))  -  (“bxay+  «y“  \*z)  "  n|Vf 


h  “  3h/3x  and  h  ■  3h/3z 
x  z 


(A.  9) 


(A. 10) 


The  conditional  densit)  can  be  expressed  as 

p(u|n)  ■  P2(n^,n^|ri>6(U-l)+[l-P2(n^,n^|n)36(U)  (A. 11) 

in  which  <S(U)  is  the  Dirac  delta  function  and  P2(iif,ni|n)  is  the 

probability  that  a  point  on  the  rough  surface  will  be  both  illuminated 

by  the  source  and  visible  at  the  observation  point,  given  the  value  of 

the  gradient  n(h  ,h  ) . 

X  z 


On  assuming  that 


k2<h2>  -  k2of  «  1,  n*a„  -  1,  P-Cn^n^n) 


(A.  12) 


the  full  wave  solutions  for  the  scattering  cross  section  (Bahar,  1981b) 

PQ 

can  be  shown  to  reduce  to  the  perturbation  solution  <0  >  (Rice,  1951; 


Barrick,  1970) 

<oPQ>  -  Air  k*W(v  ,v  )|cos6^  cos0*  aP^|2 


O  '  X'  z 


(A. 13) 


in  which  61  and  Qf  are  the  angles  of  incidence  and  scatter  with  respect 
oo 

to  the  reference  plane  (y*0) , 


W(VV2>  -  -j  J  exp(ivxxd+  ivzzd)  h(x,z)h,(x'(x')  dxddzd  (4.14) 

Is  ch«  spectral  density  of  the  rough  surface  height  (Barrick,  1970; 

Ishlmaru,  1978).  It  is  assumed,  In  this  work,  that  the  rough  surface 
is  statistically  homogeneous  and  isotropic,  thus,  the  surface  height 
correlation  function  <h(x,z)h* (x* ,z ')>  depends  only  on  distances  in 


the  mean  plane  where  <h>  «  0 


Furthermore, 


(x-x')a  +  (z-z')a  -  x  a  +  z.a 

x  z  a  a  z 


1 2  cos8^  cos©*  aP<^|2  ■  (cos0^+  cos0*)DP^  2 


(4.15) 


(4.16) 


Thus  to  obtain  the  perturbation  solution  the  unit  vector  n  (normal  to 
the  rough  surface)  is  replaced  by  a^  in  the  full  wave  expression  for  the 
scattering  coefficient  DPt^  (Bahar,  1981b). 

The  Physical  Optics  solution  for  the  scattering  cross  section 
<aP<^>  can  also  be  derived  directly  from  the  full  wave  solutions.  Thus 
(Bahar,  1981c)  for  P-V  or  H  and  Q-H,V  (P^Q) 
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(4.17a) 


(4.17b) 


in  which  tpla  la  the  angle  between  the  plane  of  incidence  associated  with 


the  reference  plane  (y*0)  and  the  local  plane  of  incidence  (associated 
with  the  local  tangent  plane)  evaluated  at  the  specular  points  where 
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Sf-  51 

|S£-  s1! 


■  v/v 


(4.18) 


Similarly  is  the  angle  between  the  reference  and  local  scatter  plane 

evaluated  at  the  specular  point.  In  (4.17)  the  Fresnel  reflection 
coefficients  and  P2(n,n1|n),  (4.12)  are  evaluated  at  6®  the  angle  of 
incidence  at  the  specular  points.  The  integral  $  is  given  by 

CO 

j  Cx2(vy,-vy)  -  lx(vy)|2]exp[ivxxd+ivzzd]dxdd2d 

—CO 

in  which  x(vy)  is  the  characteristic  function  (4.5)  and  X2(vy,-jvy)  is  the 

joint  characteristic  function  (transform  of  the  height  probability  density  p(h,h')) 


(4.19) 


w 

X2(vy.-vy)  -  |  exp[ivy(h-h,)]p(h,h,)dhdh' 


(4.20) 


The  full  wave  approach  has  also  been  applied  to  rough  surfaces  with 

two  or  more  roughness  scales.  Thus,  if  it  is  assumed  that  the  composite 

surface  can  be  represented  by  the  superposition  of  H  uncorrelated 

surface  heights  h  (n-1,2. . .N)  and  that  the  principal  feature  that  dis- 
n 

tinguishes  the  individual  surface  heights  from  each  other  is  their 


respective  correlation  distances  £n  such  that 


ll  K<  *12  <<  *2  <<  *23  <<  *3 


*N-1  <<:  *H-1,N  <<  *N 


-1  N-l  ,N 

the  scattering  cross  section  can  be  expressed  as  a  weighted  sum  of  the 
individual  scattering  cross  sections  (Bahar ,  1981c) 


(4.21) 
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(4.22) 


in  which  <aP<^>_  is  given  by  (4.3)  with  h,h'  and  X  replaced  by  hQ,h^  X° 


(the  characteristic  function  for  surface  hn>.  The  weighting  function  is  given 


by  the  product 


-jtt 


,m-l  1 2 
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in  which  x°  =  1-  Since  ) X° I  f  1  the  weighting  function  is  in  general  less 


than  unity. 

For  a  two  scale  composite  surface  such  that 

2  2  2  2 

k  <h,>  «  1,  k  <hl>  »  1 

Ol  0  2 

5  *  ‘y  '  *1  «  *12  «  *2 

the  scattering  cross  section  is  given  by  (Bahar,  1981c) 

<0PQ>  a  <0PQ>1  +  <0PQ>2 
PO 

In  (4.25)  <a  is  given  by  the  perturbation  solution  (4.13)  which  accounts 

PO 

for  Bragg  scatter.  Furthermore,  <o  >2  is  given  by  the  Physical  Optics 
solution  (4.17).  For  a  surface  with  a  Gaussian  distribution  in  (4.17)  is 
given  by 


(4.24a) 

(4.24b) 

(4.25) 


tf' 


2v 


Ik  (cos0^+  cos0^)J 
v  o  O  O' 


p(h  ,h  ) 
r  xsp  zsp 


(4.26) 


in  which  the  joint  probability  density  function  for  the  slopes  at  the 

specular  points  is  p(h  ,h  )  (Bahar,  1981c).  Thus  the  full  waves  solution 

xsp  zsp 

(4.3)  accounts  for  both  Bragg  scatter  as  well  as  specular -point  scatter. 

In  Section  4.4  the  full  wave  approach  is  applied  to  composite  rough 
surfaces  with  moderate  slopes  (h^  =1,  h z  =1). 

4.4  Scattering  Cross  Sections  for  Composite  Rough  Surfaces  with  Moderate  Slopes 
Consider  here  a  composite  surface  with  two  statistically  independent 
roughness  scales  such  that  the  radius  vector  to  the  surface  is 

r  -  ?j,(x  ,  hy.z)  +  r^h^  (4.27a) 

in  vhich 


2  2 

VV  »  1  ’  “Tx  “  ^  C  hx  ’  “Tz  “  bp*  =  hz 


(4.27b) 


fR<V  ■  \  ” 


(4.27c) 
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and 


V.2q<^>  «  1. 


dx 


1^1  «  i. 


^1  "  1^,1  «  1 


In  view  of  (4.27)  and  (4.9)  the  normal  to  the  surface  is 


n 


(4.27d) 


(4.28) 


However,  unlike  the  composite  surface  considered  in  Section  4.3,  the 
variances  of  the  slopes  are  not  assumed  to  be  small.  In  addition,  it  is 
assumed  that  the  radii  of  curvature  of  the  surface  hp  (and  therefore  h) 
is  large  compared  to  the  electromagnetic  wavelength  X.  Thus  the  full 
wave  expression  for  the  scattering  cross  section  for  the  composite  sur¬ 
face  Is  (4.3) 


-  w-  \  ^  xV 2' 

y  ■*  «-  3  n’a  -* 


•expfiv  (x-x’)  +  iv  (z-z’)]dxdz  dx’dz' 


(4.29) 


F  R  - 

In  (4.29)  x  and  x  are  the  characteristic  functions  for  surfaces  hp  and  h^ 

R  — 

respectively  and  x2  ls  the  Joint  characteristic  function  for  the  surface  h^. 

In  view  of  the  above  assumptions  for  the  small  scale  surface* 


R,  \ 
X  (v-) 


and 


<exp[iv_hR]>  *  1  -  ~  vl  <h^>  «  1,  v_  -  n-v 
<exp  iv-y(V^)>  *  1  -  \  +  <^2>)  + 

■  ixRl2  -  i  <w> 


|XF(vy)|2  «  1 

Following  the  analytical  (physical  Optics)  procedure*  presented  by 
Sancer  (1969)  for  the  large  scale  surface  hp 


(4.30a) 


(4.30b) 

(4.31) 


*The  surface  h^  which  consists  of  the  small  scale  spectral  components  cannot 
be  treated  as  the  large  scale  surface  hp  since  it  does  not  satisfy  the  radii 
of  curvature  criteria. 
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;PQ 


PPQ(r)U(r) 


n‘*. 


and 


exp  iv  (hp-h^)  exp  iv  (hxxd  +  h^) 


y'  r  r'  *  _  y- 
Thus  (4.29)  can  be  expressed  as  follows: 


<C,PQ>  ‘  nr  1  *1^  “pti»  (h  *  Hi  *  ) ] ( lx* 1 2+  V?  k^)> 
y  ln*a  1  J 


(4.32) 

(4.33) 


in  which 


<a 


PQ> 


f  < 

*  J  5.; 


•exp[ivxxd+iv2zd]dxdd2ddxdz 
1 

2 


PQ  .  PQ 
■  <0  >  +  <o  > 


(4.34) 


expC ivy (hxxd+hz  ) ] | xR 1 2> 


•exptiv  x  +  iv  z  ]dx.dz .  -  |xR|2<oPQ> 


(4.35) 


x"d  zddd 

Thus  <0Pq>o  is  the  Physical  Optics  solution  <0Pq>,  (4.17)  multiplied  by 
|X  (v*ng)|  <  1,  the  characteristic  function  for  the  small  scale  rough 

surface  evaluated  at  v*ng  the  projection  of  v  on  the  normal  at  the  specular 


point.  Furthermore, 

<aPq>  -  jLf 

1  it  J 


DPQUv_ 


n’a 


eXP[i(Vx'>Vyhx)Xd+i(VZ+VyhZ)2d:|dXddZd>  (4>36) 
PQ 

Thus  <a  reduces  to  the  perturbation  solution  (4.13)  when  the  surface 

slopes  are  small  (Bahar,  1981b).  For  the  case  considered  in  this  section  the 

distribution  of  the  large  scale  surface  slopes  has  an  effect  on  the  component 

PQ 

of  the  scattering  cross  section  due  to  the  small  scale  roughness  <o  >^. 

Note  that  the  small  scale  roughness  h^  is  a  perturbation  about  the  large 
scale  roughness  and  h^  is  the  perpendicular  distance  from  the  unperturbed 
(filtered)  surface  y  -  hpOt.hy.z)  to  the  perturbed  surface  (4.27) 

(Burrows,  1973).  Thus  the  expectation  of  hR(x,z)h^(x’ ,z ’)  is  a  func¬ 
tion  of  distances  measured  along  the  plane  tangent  to  the 
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1 


surface  y  -  hp(x,z)  -  0  (Wright,  1966;  Valenxu«la,1968)  and  not  dis¬ 
tances  measured  along  the  mean  plane  lrnl  (4 .15) .  To  facilitate  the 
evaluation  of  (4.36)  the  integrand  is  expressed  in  terms  of  the  unit 
vectors  n^  and  n^  and  variables  x,y,z  of  the  local  coordinate  system 

associated  with  the  local  tangent  plane  (Bahar,  1981a;  Bahar,  1981b). 
The  unit  vectors  n3>n2,n3  can  be  chosen  as  follows  (see  Fig.  4.2) 

n.  -  (n  x  a3)/|n  x  a3|,  ii2  -  n,  n3  ■  n3  x  n  (4 


(4.37a) 


in  which 


a,  ■  a  ,  a_  ■  a  ,  a,  *  a 

1  x  2  y  3  z 


(4.37b) 


-1  -  -2  -  -3  - 

and  the  relationship  between  the  local  coordinate  system  X(x  -x,x  -y,x  -z) 

12  3 

and  the  reference  coordinate  system  X(x  “x,x  «y,x  «z)  is 


in  which 


Thus  in  (4.36) 


x1-  ■  a^x^  (sum  on  j-1,2,3) 


aij  “  V*j  i’j“1*2*3 


(v  +v  h  )x,+(v  4v  h  )z.  ~  v-  X  +  V-  z 
X  y  X  d  z  y  Z  u  XU  ZQ 


in  which  v_,  v_,  v_  are  the  components  of  v  (4.4)  in  the  local 
x  y  “ 


coordinate  system 


v_  -  vnr  -  v-n2,  vf  -  v*n3 


Furthermore  (see  Fig.  4.3) 


xd  -  x-x '  ,  ld  -  z-z ' 


(4.38a) 


(4.38b) 


(4.39a) 


(4.39b) 


(4.40a) 


■V'd 


'  dh  did 


Therefore  (4.36)  can  be  expressed  as 
pn  . 


/Q>  -  ~  <j— -Z-  «pI>sV  iv**d]dxd  dxc 

>  n*a 


(4.40b) 


(4.41) 


r*  =  DISTANCE  MEASURED  ALONG  LARGE  SCALE  SURFACE 

3  Rough  surface  decomposed  into  large  and  small  scale  surfaces. 
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Furthermore  assuming  that  the  statistics  of  the  small  scale  surface 

are  independent  of  the  slopes  of  the  large  scale  surface  h_*  (Valen- 

F 

zuela,  1968,  Brown,  1978)  (see  Fig.  4.3) 


JL 

2 


<  1^  (x , z ) h^ (x ' , z ' ) >exp ( iv_xd+iv^zd ) dxd  d zd 


W(v_,v_) 
x  z 


(4.42) 


IT  ' 

in  which  the  dependent  variables  of  the  spectral  density  of  the  small 

scale  roughness  v_  and  v_  are  the  orthogonal  components  of  v  in  the  local 
x  2 

tangent  plane.  Thus  W(v_,v^)  depends  not  only  on  the  direction  of  the 
incident  and  scattered  waves  but  also  on  the  normal,  n,  to  the  rough 
surface.  Substitute  (4.42)  into  (4.41)  to  get 


|DP(W'2 


<aP<^>.  *  TTk^f  - W(v_,v_)p(n)p(U  |n)dh  dh  dlJ 

1  of—  —  x  z  xz 

'  n*  a 


(4.43) 


y 

and  substitute  (4.11)  for  P(U|i5)  unto  (4.43)  and  integrate  with  respect  to  U 
to  get 


<0 


PQ>_ 


TTk 


,-f  riirs 


^ —  W(v-,v-)P7(ni ,nA|n)p(n)dh  dh 


n*a 


Note  that  the  conditional  probability  P2(nf,ni|n)  for  arbitrary  slope 

n (h  ,h  )  can  be  expressed  in  terms  of  the  conditional  probability 
x  z 

P  (n^,ni|n  )  for  the  specular  points  (n  -*■  n  )  as  follows  (Smith,  1967; 
2  S  S 

Brown,  1980) 


(4.44) 


P2(nf,ni|n)  •  S(n*n^)S(-n-ni)P2(nf ,n* |ns)  (4.45) 

in  which  S(a)  is  the  unit  step  function.  The  arguments  of  S.-n’n1  and  n>nf 
vanish  at  points  of  the  rough  surface  where  the  incident  and  scatter  wave 
normals,  n"1  and  n^  are  tangent  to  the  surface. 

In  order  to  interpret  the  results  (4.44)  assume  that  surface 
h_  is  normal  to  the  constant  vector  nQ 


"This  assumption  of  statistical  independence  of  the  small  scale  R 

R 

and  slopes  restricts  the  class  of  rough  surfaces  assumed  in  this  work. 
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n  -*•  n  ■  n  a  +n  a  +  n  a 
o  ox  x  oy  y  oz  z 

The  equation  for  this  surface  is 

n-r-n  x+n  y+n  z  ■  0 
o  ox  oy  '  oz 

thus  for  y  -  hF(x,a) 

h_  «*  -(n  x  +  n  z)/n 
F  ox  oz  oy 

Furthermore 


(4.46a) 


(4.46b) 


(4.46c) 


p(n)  »  6(h  +  n  In  )6(h  +  n  In  )  •  6(n-n  )  (4.47) 

e  x  ox1  oy  z  oz 1  oy  o 

Assuming  that  n*n  >0  and  -n-n*  >  0  (the  plane  n  *r  ■  0  is  visible  and 

o 

PO 

illuminated),  the  normalized  scattering  cross  section  <a  (4.45)  is 
given 


<a 


pQ>. 


TTk 


W<V5c*Vg)  l^^yl2 


n  *a 
o  y 


(4.48) 


n=*n 


Thus  the  scattering  cross  section  (4.48)  for  the  inclined  plane  (4.46)  is 
in  complete  agreement  with  the  perturbation  result  (4.13)  for  arbitrary 
nQ.  It  is  equal  to  it  for  nQ  »  a^  (Rice,  1951;  Barrick,  1970).  The  term 
appearing  in  the  denominator  of  (4.48)  is  a  consequence  of 
normalization  (4.2)  (since  A^  is  the  projection  of  the  surface  A  on  the 
reference  plane) . 

The  integrand  in  (4.44)  is  therefore  the  contribution  to  the  scattering 
PQ 

cross  section  <a  from  "patches"  of  the  rough  surface  that  are  normal 
to  the  vector  n  (Wright,  1968;  Valenzuela,  1968). 

In  his  review  of  theoretical  treatments  of  the  scattering  of 
electromagnetic  waves  by  rough  surfaces,  Valenzuela  (1968)  states  that 
scattering  models  have  been  developed  "that  make  it  possible  to  apply 
available  theories  to  surfaces  that  cannot  be  treated  purely  by 
perturbation  and  Physical  Optics.... in  these  models  composite  surface 
approximate  the  rough  surface  by  'patches’  of  slightly  rough  surfaces 
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thac  ride  the  large  waves.... the  total  scattered  power  from  this  composite 
surface  is  obtained  by  averaging  over  the  distribution  of  slopes  of  the 
large  waves."  Valenzuela  notes  however,  that  "these  composite  models 
are  mostly  based  on  physical  considerations  and  are  able  to  explain 
features  in  radar  cross-section  data  from  the  sea  that  no  theory  can." 

It  is  shown  in  this  section  that  using  the  full  wave  approach,  one  can 
provide  a  self-consistent  theoretical  basis  for  scattering  by  composite 
surfaces  without  using  a  combination  of  Physical  Optics  and  perturbation 
theories. 

To  compare  this  work  with  earlier  analysis  of  the  backscatter  cross 
section  (Valenzuela,  1968;  Valenzuela,  Laing  and  Daley,  1971)  in  which 
the  integral  (4.44)  is  expressed  in  terms  of  tilt  angles  <p  and  6  "in" 
and  "perpendicular"  to  the  reference  plane  of  incidence  (normal  to 
n^  x  Sy)  rather  than  in  terms  of  the  slopes  h^  and  h^ ,  the  normal  n 
to  the  rough  surface  is  expressed  as  follows  for  n*  in  the  (x,y)  plane 

n  ®  sinijj  cos6  a  +  cos4<  cos6  a  +  sin6  a 
x  y  z 

-  (-  h  a  +  a  -  h  a  )/(l  +  h2  +  h2)**  (4.49) 

x  x  y  z  z  x  z 

In  the  work  by  Valenzuela,  et  al.  (1971),  shadowing  had  been  ignored, 
however  the  arguments  of  W(v_,v_)  in  their  work  is  in  agreement  with  (4.44). 

The  result  (4.44)  is  also  in  agreement  with  the  expression  derived 
recently  by  Brown  (1978)  (as  corrected  in  Brown,  1980)  provided  that  the 
slope  of  the  large  scale  rough  surface  is  very  small.  The  differences 
between  these  two  solutions  are  primarily  due  to  the  assumption  in  Brown's 
work  that  h^  (£g  in  Brown's  notation)  is  the  distance  from  the  unperturbed 
(filtered)  surface  to  the  perturbed  surface,  measured  along  a  line  per¬ 
pendicular  to  the  reference  plane  (y-0)  rather  than  perpendicular  to  the 
unperturbed  surface.  This  is  contrary  to  the  original  analysis  of  Burrows 
(1973)  upon  which  Brown's  analysis  (1978)  is  based.  Furthermore  in  this  work, 
the  expectation  of  h^h^  is  considered  to  be  dependent  on  distances  r^  measured 
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along  the  plane  tangent  to  the  surface  y-h^x.z),  (i.e.,  the  surface 

normal  to  n)  where  r^  *  xd  ni  +  zd  n3*  *n8tea<*  distances  measured 

along  the  mean  (reference)  plane  r  ■  x,  a  +  z,  a  (4.15)  as  assumed 

to  a  x  <2  z 

by  Brown  (see  Fig. 4. 3) .This  feature  makes  the  full  wave  solutions 
invariant  to  coordinate  transformations  (4.48).  In  his  analysis. 

Brown  (1978,  1980)  accounts  for  shadowing  and  uses  both  Physical  Optics 
and  Burrows  (1973)  perturbation  theories  to  derive  his  results  in 
terms  of  a  two  dimensional  convolution  of  transforms. 

4.5  Concluding  Remarks 

Composite  models  of  random  rough  surfaces  with  moderate  mean  square 
slopes  are  analyzed  using  the  full  wave  approach.  The  full  wave  approach 
is  shown  to  account  for  specular  point  scattering  and  Bragg  scattering 
in  a  self ^consistent  manner  without  introducing  Physical  Optics  and 
perturbation  theories.  Thus  using  the  full  wave  approach,  it  is  not 
necessary  to  specify  the  wavenumber  k^  where  the  rough  surface  spectral 
splitting  is  assumed  to  occur  (Brown,  1978,  1980).  However,  decomposi¬ 
tion  of  the  surface  into  a  large  scale  surface  h^  and  small  scale 
surface  h^  enhances  one's  physical  insight  and  facilitates  the  numerical 
evaluation  of  the  scattering  cross  sections.  The  wavenumber  kj  is  used 
to  distinguish  between  the  large  scale  rough  surface,  h^,  and  the  small 
scale  rough  surface,  hg.  The  solutions  derived  in  this  paper  are 
compared  with  solutions  derived  earlier  "mostly  based  on  physical 
considerations"  (Valenzuela,  1968).  They  are  also  shown  to  be  in 
agreement  with  Brown's  solution  provided  that  the  mean  square  slope 
of  the  rough  surface  is  very  small  (Brown,  1978,  1980).  Brown  uses  a 
combination  of  Physical  Optics  (Beckmann,  1968)  and  Burrows  (1967) 
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perturbation  theories  to  develop  his  solution.  The  differences  between 
Brown's  solution  and  the  full  wave  solution  (4.44)  are  primarily  due  to  the 
assumption  in  Brown's  work  that  ^  is  the  distance,  from  the  unperturbed 
surface  to  the  perturbed  surface,  measured  along  a  line  perpendicular  to  the 
reference  plane  (y-0)  rather  than  perpendicular  to  the  unperturbed  (filtered) 
surface.  This  is  contrary  to  the  original  analysis  of  Burrows  (1973)  upon 
which  Brown's  (1978)  analysis  is  based.  Furthermore,  Brown  assumes  that  the 
autocorrelation  function  for  the  small  scale  surface  height  depends  upon 
distances  measured  in  the  mean  plane  rather  than  distances  along  the  plane 
tangent  to  the  large  scale  surface  as  assumed  here  and  implicitly  by 
Valenzuela  (1968). 

The  contribution  to  the  scattering  cross  section  due  to  Bragg 
scattering  (associated  with  the  small  scale  surface)  is  given  by  an 
integral  over  the  slopes  of  the  rough  surface.  This  contribution  is 
proportional  to  and  depends  on  the  polarization  (4.44).  It  is  shown 
to  be  in  complete  agreement  with  perturbation  theory  (Rice,  1951; 

Barrick,  1970;  Barrick  and  Peake,  1968).  The  integrand  in  this  expression 
is  proportional  to  the  small  scale  surface  height  spectral  density 
W(v^,v_)  whose  arguments  are  the  components  of  v  •  kQ(n^-n*)  in  the  local 
tangent  plane.  Thus  W(v_,v_)  is  also  a  function  of  slopes.  The  contribu¬ 
tion  to  the  scattering  cross  section  due  to  specular  point  scattering 
(associated  with  the  large  scale  surface)  is  independent  of  kQ  and 
polarization. 

Using  the  full  wave  approach  the  specification  of  kj  (where  spectral 

splitting  is  assumed  to  occur  is  not  restricted  by  the  perturbation 
2  2 

condition,  (5  -  4ko<hs>  «  1.  Thus  it  can  be  applied  to  more  general 

models  of  rough  surfaces  for  which  a  combination  of  solutions  based  on 

perturbation  and  physical  Optics  are  not  suitable.  The  specifications  of 

PQ 

k.  and  the  sensitivity  of  the  cross  sections  <0  to  variations  in 
d 

k.  are  subjects  of  future  investigations, 
d 


v-  *'•'  • 
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5.0  SUMMARY  OF  RESEARCH  TO  BE  CONDUCTED  DURING  SECOND  YEAR  OF  CONTRACT 
5.1  Analysis 

In  order  Co  account  for  both  specular  point  scattering  as  well  as 

Bragg  scattering,  a  two  scale— composite  model — of  the  rough  surface 

is  generally  used.  The  specular  point  scattering  is  accounted  for 

through  the  use  of  Physical  Optics  theory,  while  perturbation  theory 

is  used  to  account  for  Bragg  scattering.  The  principal  difficulty 

with  this  perturbed-Physical  Optics  approach  to  analyze  composite  models 

of  rough  surfaces  lies  in  the  specification  of  the  wavenumber  kj  where 

the  surface  height  spectral  splitting  is  assumed  to  occur.  On  one  hand  kj 

must  be  sufficiently  small  such  that  the  radii  of  curvature  of  the  large 

scale  (filtered)  surface  is  large  enough  to  justify  the  application  of 

Physical  Optics  theory  to  the  filtered  surface.  On  the  other  hand  k. 

a 

must  be  sufficiently  large  such  that  the  mean  square  height  of  the  small 
scale  surface  is  small  enough  to  justify  the  application  of  perturba¬ 
tion  theory  to  this  small  scale  surface.  In  general  these  two  restrictions 
on  the  specification  of  k^  may  not  be  satisfied  simultaneously.  It  is 
found  using  this  perturbed-Physical  Optics  approach,  that  the  computed 
value  of  the  rough  surface  scattering  cross  section  will  critically  depend 
upon  the  specified  value  of  k^.  Using  the  Full  Wave  approach  (which 
accounts  for  both  specular  point  scattering  as  well  as  Bragg  scattering 
in  a  self-consistent  manner),  the  specification  of  the  wavenumber  k. 

U 

(at  which  spectral  splitting  is  assumed  to  occur),  is  not  restricted 
by  the  conflicting  considerations  of  Physical  Optics  and  perturbation 
theories.  Rather  in  this  case  k^  is  determined  only  by  analytical  and 
computational  considerations.  Thus,  if  kj  is  chosen  judiciously  the 


computed  value  for  the  resulting  scattering  cross  section  should  not  depend 


upon  the  specific  value  of  k,  chosen  for  these  computations.  Instead  k,  is 

a  d 

determined  solely  by  the  need  to  minimize  the  numerical  computations. 

This  aspect  of  the  problem  is  being  investigated  in  detail. 

During  the  next  year,  the  Full  Wave  approach  will  also  be  applied  to 
non-Gaussian  rough  surfaces.  Thus  the  following  topics  will  be  investi¬ 
gated  in  detail: 

a.  Statistical  description  of  classes  of  non-Gaussian  rough  surfaces. 

b.  Shadow  functions  for  non-Gaussian  rough  surfaces. 

c.  Like  and  cross  polarized  scattering  cross  sections  for  non- 
Gaussian  rough  surfaces. 

Rough  surfaces  for  which  decorrelation  implies  statistical  independence 
are  of  particular  interest  in  our  investigations. 

5.2  Computer  Program 

The.  computer  programs  used  to  evaluate  the  scattering  cross  sections 

for  composite  rough  surfaces  are  being  modified  to  reflect  the  recent 

advances  made  in  the  analytical  approach  to  this  problem.  The  effect  of 

the  choice  of  k^  on  the  numerical  computations  Is  being  investigated. 

From  the  preliminary  results  it  is  shown  that  the  numerical  value  for  the 

scattering  cross  sections  as  determined  by  the  recently  modified  computer 

programs,  are  not  dependent  on  k.  and  that  the  most  suitable  choice  for 

a 

2  -2 

k,  is  such  that  4k  <h>  =  1  (k  is  the  electromagnetic  wavenumber  and 
a  o  K  o 

-2 

<h£>  is  the  mean  square  of  the  small  scale  surface  height).  More  work 
needs  to  be  done  to  complete  this  phase  of  the  investigation. 


